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1 . INTRODUCTION 



This report represents the second phase of a long-range project 
to study self-adaptive control principles and apply them to the 
control of aircraft. The primary purpose of this phase is to develop 
a suitable mathematical model for the study of aircraft dynamics 
including all degrees of freedom, and the secondary purpose is to 
test the feasibility of using active compensators in the autopilot 
as a basis for self-adaptive schemes of control. 

It was decided to restrict the model to that of a rigid airframe, 
and it was also decided to use a block diagram model because of 
well established correlations between the block diagram and a 
determinantal array on one hand, and between the block diagram 
and analog simulation on the other hand. The initial block diagram 
model IS to a linear model for two reasons: first, it is necessary 
that the model be reducible to simpler, well-known linear block 
diagram models when the proper restrictions are imposed; second, 
analysis of a linear model is readily obtained and provides a basis 
for estimating the importance of nonlinear effects as they are 
introduced. It is anticipated that nonlinear effects can be introduced 
by the injection of properly derived signals, and it is expected that 
this modification of the model will leE^d to simulation studies since 
the complications caused by the nonlinearities may make an analytic 
treatment impossible. 
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Another reason for choosing a block diagram model for the 
airframe is the ease with which an autopilot configuration can be 
added to the model* It is expected that the self-adaptive compen- 
sation devices will be incorporated in the autopilot, though it is 
believed that the block diagram model may lead to profitable 
analysis of the aerodynamic characteristics of the airframe, 
with resultant indications of ways to improve its design. 
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2, DEVELOPMENT OF THE STANDARD BLOCK DIAGRAM 
AND THE STANDARD DETERMINANT. 

2.1 Piirpose 

it IS the puppaae e£ this seetien to develop the mathematical 
background leading to the derivation of the block diagram model. 
It will be shown that there is a unique relationship between a 
properly arranged block diagram and a determinantal array such 
that either may be wi'itten down by inspection of others. The 
theoretical development starts with the block diagram and leads 
to the mathematical representation; when applied to the airframe 
one starts with the simultaneous equations and arrives at the 
block diagram. 
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2 . 2 The Standard Block Diagram 

The theory of simple, single loop control systems is reason- 
ably well developed,, and the block diagrams of such systems are 
n©l On ©the? hand the thee?y o£ multiple -loop, 

multiple -input, multiple -output and multicoupled control systems 
is not well developed because of the complexity of the problems 
involved. The block diagrams of such systems are naturally 
complicated and sometimes difficult to interpret. The obvious 
approach to the analysis of such a system is the formulation and 
solution of a set of simultaneous differential equations for the 
system. Perhaps the easiest way to formulate the equations is 
from a block diagram of the system. In using this approach 
certain advantages accrue from standardizing the block diagram, 
and certain other advantages accrue from generalizing the diagram 
By standardizing is meant specifying the arrangement of the sum- 
ming points, blocks and pick- off points, and specifying a notation 
for the transfer functions. By generalizing is meant extending 
the standardization rules so that the maximum number of blocks 
is clearly specified and organizing the mathematical approach so 
that the presence or absence of any block is readily evident in the 

V 

equations. It is clear that few systems will contain, the maximum 
number of blocks, so that the ability to recognize the presence or 
absence of blocks from the equation is in itself an advantage. It 
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is also clear that block diagrams set up from physical analysis 
must be manipulated to gef them into the standard form. 

Any system contains* a finite number of summation points, 
earn© ©f th§§# in a nife.in tt‘angjrni§&ion pfttH in th§ 
that the signal output from the summer is fed into a block which 
produces a signal output closer to a system output. If all other 
summing points are ignored these main summing points are 
?i.-’/ttered in sequence, a, b, c, d . . n, and the standard block 
diagram is drawn by arranging these summing points in a straight 
line from left to right, with intermediate blocks 'as required. It 
is obvious that the numbering sequence is not completely arbitrary. 
In general summer ”a” L^eds a block whose output feeds summer 
”b" , etc. , though there may be breaks in the signal flow as will 
be shown later. Also the output of the block following the nth 
summer, is normally a system output though additional system 
outputs may be derived at preceding points. 

Each summer may receive a number of signals, but delivers 
only one output. At each summer there is the possibility that 
one of the signals received is an input command from an external 
source or a disturbance signal such as noise. No distinction is 

made here. It is convenient to designate such inputs (at their 

' '■ ' **' 

sources) by numbers in sequence 1, 2, 3, etc. It is also corirven- 
ient to relate the input sequence with the summer sequence sp 
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< 

that input 1 feeds summer "a", input 2 feeds summer ''V\ ^tc . 

, I 

If any summer does not have an external input signal, the 
corresponding number for the input signal is not used. In writing 

the aquationi the signals from ths input stations ars designated 

by Roman numerals- 



The, .transfer functions designated in the blocks shall, -in 
general, be designated by double subscript notation. To imple- 
ment this the signal at the output of a summer is designated by 
the same lower case letter which designates” the summer, and 
the signal output of the block which receives this signal is desig- 
nated by the same letter, but using the upper case form. This 
has the advantage of clearly indicating all important signal sources 
and the disadvantage of permitting two symbols which represent 
the same transfer function, as will be shown. The double sub- 
script notation shall consist of two letters or a combination of a 
number and a letter. The first subscript designates the input to 
the block and the second subscript designates the, summer into 
which the output of the block is fed. 

Thus , 



G, = .transfer function of a block between input 1 and summer "a 
la 

G_ = transfer function of a block between input 1 and summer ”b 
' lb ^ 

G . = G , = transfer function of a block between summer "a" 
aA ab 

and summer "b" 
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G = transfer function of a feedback block between signal 
^^a • 

A arid summer ’’a" ' . 

= transfer function of a feedback block between signal 
Ba ^ 

B and summer ”a” 

. G^^= transfer function of a feedback block between signal ’ 

C and summer "a” 

Many other interconnections are obviously possible and are 
readily handled by the double subscript notation. 

It is common terminology to speak of systems in terms of the 
numbers of loops involved. From the viewpoint of mathematical 
analysis the number of loops is not a convenient description of 
the complexity of the system. A better ‘Criterion is the number 
of summing points or ’’nodes," and this terminology will.be used 
' in this text . 

.The development or manipulation of the standard block dia- 
gram is most readily shown by a’ series of illustrations. Figure 2. la 
shows the simplest case of a single node system. For many 

practical systems G = G = 1.0, All other variations of the 

1 a ir\.a 

single node' system may be reduced to this form if desired 
Figure 2.1b shows a single node system with additional blocks and 
summers representing the insertion of tachometer feedback 
and input velocity feedforward. Note that the simple rules for 
block diagram manipulation permit reduction of Figure 2,1b to 
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a) The Si ai. -lard Diagram 




b) A block diagram which readily manipulates to the ctuadard form 
Fig. 2-1 STANDARD BLOCK DIAGRAM FOR SINGLE NODE SYSTEMS 




Fig. 2-2 A TWO NODE SYSTEM DERIVED FROM FIG. 2- lb. 



the form of Figure 2. la. Of course, it is not always desirable 

to carry out^ such a reduction. If the particular problem involves 

determination of suitable parameters for G , a slight manipula- 

B b 

tion rasults in the block diagram of Figure 2.2, which is a two 

node system. In this .figure G is a new transfer function obtained 

’* • 13 . 

by manipulating Figure 2.1b, and G^^= 1.0. 

The block diagram of Figure 2.2 is in the standard form, but 
is a special case since it obviously does not contain all the signals 
and blocks possible for a two node system. A generalized block 
diagram for a two node system is shown in Figure 2.3. Note that 
each summer receives a number of signals, and it is necessary 
to designate whether the signal adds (+) or subtracts (-). The 
conf.igur.ation* of Figure 2. 3 is easily expanded to the block diagram 
of a n-node system as shown in Figure 2,4. 

r 

I 

The preceding figures have illustrated the multiple node block 
diagram for cases with only a single output but many inputs. In 
practice any number of outputs may exist in combination with one 
or many inputs. The basic procedure for putting the block diagram 
is standard form is unchanged, but the results of the manipulations 
sometimes take a peculiar configuration. A number of block dia- 
grams are therefore presented to illustrate some of the possible 
configurations . 

’ Figure 2.6 shows a system with one input and two outputs. 

r]' ■ ^ ^ ■ • 

The summers are lettered and arranged in sequence as usual. 
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II 




■ Fig. 2.3 GENERALIZED BLOCK DIAGRAM FOR A TWO NODE SYSTEM 



10a 



4 




Fig. 2.4 STANDARD BLOCK DIAGRAM OF AN N-NODE SYSTEM 






p 

p. 
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Fig. 2.5 ONE INPUT; TWO OUTPUTS 



There is a break, or cut, in the main path since the output of 

G . is not fed into node b.' Figure 2.6 shows a system with 
3./V 

t 

two inputs, two outputs, and a teduced, number of paths coupling 

iHitnedesi ' 

i > 

2 . 3 Mathematical Analysis Based on the Standard Block Diagram 
The obvious approach to system analysis is to derive a set 

4 

of simultaneous equations for the system. This is readily accom- 
plished by writing an equation for the signal summation at each 
summer or node. If this is done for the n-node standard block 
diagram of Figure 2,4 (with proper attention to* subscripts in 
arranging the terms) a determinant can be obtained which has a 
sort of symmetry. In addition, as will be shown, each term in 
the determinant has a specific relationship to a definite path in 
the system so that a reasonable amount of physical correlation 
exists. 

Referring to the n-node block diagram, note that at each 
node there is only one output signal, but many signals being fed 
in. Of the signals being fed in some are from input stations, 
some are feedforward signals, and some are feedback signals. 
From the nature of a summing device it follows that at each node 
the output signal must be the algebraic sum of all input signals. 
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TT 



2 



1 




output 

#1 



Fig. 2-6 TY/O INPUTS; TWO OUTPUTS 
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In feedback ^control systems inverse or negative feedback is 
normally used, so it is, convenient to treat feedback signals as 
subtrftGtiVCi Input and feedforward lignals may either ad or 
subtract, so it is convenient to consider them additive in the 



generalized treatment. 

At node "a" the output of the node is: 






(2-1). 



(G G^)d-...-(G G )n+G, I + G^ II 
de Da no Na la 2a 



+ m + G. IV + . . . G 



3a 



4a 



na 



H • 



This rearranges to 



{ 1 + G , G^ ) a +(G, G^ ) b + (G , G^ ) c + 
ab Aa be Ba cd Ca 



+ (G G_ ) d + ... (G G .) n = 
' de Da 'no Na 



m = n ( f/f ) 



m = 1 (I) m; 



(G )m (2-2) 



At node "b" the arrangment of signals is similar to that at 

node "a", but the first feedforward signal appears. It is the 

signal fed forward from node "a" through G , . The equation at 

a [/ 

node "b" may be written as: 



» = l=ab> “ - <°bc °Bb> - <°Cb - Sb> " 

- <°„o °Nb> “ * "=lb ' ^ °2b “ " Sb “ " ■••^°nbH 
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Keeping the signals in normal sequence this rearranges to; 



^ <°cd °Cb> " °de °Db> ^ " 



...+(G G )n = 
no I'lb 



m ) 

m = 1 (I ) 



(2-4) 



<°mb> 



In’ like manner the equilibrium equation at node "c" is: 



- = <°ac> " " <°bc> » - <°cd °Cc> <°dc °Dc> 



..-(G G^ In+G, I + G, II G, III + . . . + G \A 
' no Nc' Ic 2c 3c • nc 



(2-5) 



This rearranges to ■ 



-(G ) a - (G^ ) b + (1 + G ^ G^ ) c + (G^ G^ ) d + 
' ac' ' be' ' cd Cc ' de Dc' 



m = n OA) 



(2-6) 



+(G G ) n = (G ) m 

no Nc . me 

m = 1 (I ) 









By inspecting equations 2-2, 4,; 6, the symmetry of arrange- 

ment of terms is easily determined, and the equation at any node 

may be written without the need to refer to the block diagram. It 

should also be noted from these equations that all of the input signals 

or forcing functions are located on the right hand side of each 

equation leaving- on the left hand side the terms which refer to 

components and signals within the closed loops. Thus the. left hand 

/ ‘ 

sides may be arranged to form the characteristic determinant of 
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the system. This determinant may be called and is: 



A = 



1 + G G 


+ G ~G 


0 

o 


G ^G^ 


G G 


a'A Aa 


bB Ba 


c C Ca 


dD Da 


. nN Na 


- G • 


1+G G_, 


G ^G„ . 

' cC Cb 


G G --- 


-G ^,G , 


ab 


bB Bb 


dD ab 


nN Nb 


-G 


-G, 


1 + G ^G^ 


G G_ -- 


- G G 


ac 


be 


cC Cc 


dD Dc 


nN Nc 


-G , 


— G_ 


-G , 


1 + G G_ . — — ' 


-G G , 


ad . 


bd 


cd 


dD Dd 


nN Nd 


• 


• 


* 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


-G 

an 


-°bn 


-G 

cn 




1+G 

nN N] 



(2-7) 

The arrangement of the characteristic determinant has certain 



features which can be useful: 

a) Every term in the main diagonal is of the form 1 
and each of these terms is itself the characteristic 
equation of the basic loop at the designated node. 

b) Every term above the main diagonal represents a feedback 
path consisting of forward transmission from a node through 
one block, then transmission from the output of that block 
through one feedback block to a prior node. The node at 
which this path starts is designated by the column, and the 
node to which the signal is fed is designated by the row. 

c) Every term below the main diagonal represents a feedfor- 
ward path through one block only. The starting point of 
each path is the output of the node designated by the 
column and the node to which the signal is fed is designated 
by the row. 
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In practical systems not all paths will exist. When a path 



docsnot exist the transfer function of the block which is omitted 
is zero. On the other hand if the path exists but there is a direct 

than transmission through a blocks then that 
transfer function is unity. Thus the characteristic determinant 
of any specific system may be obtained easily by arranging its 
block diagram in the standard form, writing the determinant 
of the general system for the given number of nodes, and inserting 
the transfer functions of the specific system at the proper location 
in the determinant structure. 

The characteristic determinant of any system is simply tha,t 
portion of the N-node determinant which is specified by the number 
of nodes in the given system. For example, for a one node system: 

lx (one node) = 1 + ^2-8) 

for a two no(ie system: 




G G 
bB Ba 



A (two nodes) = -G , 

ab 




(2-9) 



for a three node system: 

A (three nodes) = 1 + G G 



aA Aa bB Bb 



^cC^^Ca 



( 2 - 10 ) 



G 



ab 




- G 



- Q 



bB 




ac 
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The solution of the system equations for system perform- 
ance is obtained with the usual manipulations of the determinant. 
That is, the column representing the input functions replaces 'the 



d6iir@4 in th@ ^ni>nif€H§fi@ 

and the resulting determinant is the numerator determinant A 
is the denominator determinant, and expansion of these deter- 
minants provides the solution for the desired signals. If all 

* 

inputs exist and are activated, then the column for the input 
functions is: 




(G ) (m) 



ma 



m = n (l/\) 
m = 1(1 ) 



( 2 - 11 ) 



m = n (V^) 




(G ) (m) 
__ me 



m = 1 (I) 



m = n ( H) • 



(G ) (ni) 
mn 



m = 1 (I) 






This column replaces the nth column in the general determiniSt 
of equation 2-7 to form the numerator determinant. When any, other 
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r 

signal is to be computed Equation 2-11 replaces the corresponding 
column in the numerator determinant. 

In most practical cases only one input is activated at a time, 

80 that thi termB in Equation 2*^11 are relatively simple# In analysis 

I 

and design computations the test functions likely to be used as 
inputs are the step and the ramp, and for these the terms in 

Equation 2-11 are extremely simple. It should be noted that this 

« 

theory has been developed for linear systems, so the principle of 
superposition is valid. Thus even in the case r>f multiple inputs 
acting simultaneously the system may be analyzed for response 
to one input at a time and the results superimposed. 
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3. AN ILLUSTRATION OF THE USE OF THE DETERMINANTAL 
APPROACH. 

3. 1 Introduction 

A high gain uncompensated positioning servo is shown in 

, I 

block diagram form on Fig. 3.1. It is the purpose of this section 
to study the effect .of amplidyne voltage feedback, either positive or , 
negative, tp the tnag amp’ inpxit. 

■Values of K and have been determined as 
K = 57.5 and = .315 



3; 2 The Characteristic Determinant and Steady-State Error 

The characteristic determinant may be written down by inspec- 
tion of the block diagram Fig. 3. 1. This determinant is; 

.3 



A= 



1 



■50 



315 X 10' 
s{s-H00)^s^+10s+l00) 



, 11,500 

1 + — G 

s + 200 c 



- 11,500 



s + 200 



Expansion yields: 

A . 11,500 18.1x10^^ (3-2) 

8 + 200 c s (s +200) (s+100)'^ (s^ + 10s + 100) 



.10 



(3-1) 



The steady state error from a step input is: 



ss 



= Lim 
step s > 0 



1 - ^cd Aic 



A 



(3-3) 
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MAGAMP AMPLIDYNE MOTOR & LOAD 




a/ 



Fig. 3.1 BLOCK DIAGRAM OF A CONTROL SYSTEM 



Evaluating each term independently before combining; 



From Equation 3-2: 

Lim = 1 + 57,5G + 905 

- ^ s 

8^0 



Lim 



ac 



s ^ 0 



50 



= 2876 



1 + 57. 5G 



-57.5 



Lim G ^ = .315 
cd 



s > 0 

Substituting into Equation 3-3 yields 



ss 



step 



-t- 
= 0 -; 



.315 



(2876) 



+ 905 1 



1 + 57.5Gc + 905 

s 



8 » 0 



905 



+ 57. 5G s + 905 _ 

c s ^ 0 






. n . - . n 

L 57.5 G s + 905 . 

c s > 0 



!• - 



CO 



905 

ri57 5s K (z_) (z^) (z ) 



- c 1 



n 



(P)) (Pj) — (p„) 



s 0 



This error will be zero if there be no p = 0. 
The steady state error for a ramp input is: 



ss 



ramp 



t- 



° 



cd 



A 



-J 



s ^ 0 






(3-4) 
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Substituting values: 



E 



ss 



ramp 



= i [l - _ 

s 1 



905/i 



] 



+ 57.5 G + 905/s 

c s ^ 0 



= J_ 
s 



1 +57 G + 905/s - 905/ 
c 

1 + 57.5 G + 905/s 
c 



3 



s + 0 



1 + 57. 5G 



8 + 57.5 G s + 905 
c 



s p 0 



E , 
ss 



ramp 



K (z ) (z ) --- (z ) 

lj.57 5 ^ ^ ^ ^ 

1 +57.5 (P3) ---l?7x) 

57.5 G s + 905 
c 



s > 0 



Compensators containing terms of the form are not allowed. 

s 



Therefore, 

E 



1 + 57.5 K 



ss 



ramp 



(^l> 

^ *Pl> <P2> — 

905 



In order that the steady state error be as small as possible, the 
term 



57.5K (.,)(z^) — (.„) 
(Pj) (P^) — (P„) 



(3-5) 



This may be best obtained by having one or more of the z’s identi- 
cally zero. Then the steady state error is the minimum obtainable 
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3. 3 Summary of Compensation Schemes V/hich Appear 

Feasible and Restrictions to be Observed When Attempting 
Compensation. 

Cornpensators containing a term 1/s are prohibited because 

in ¥@du^ing iHe snvvn In ^ type aei>n 
system and producing a K =0 and an infinite error to a ramp 
input. 

In all other compensators, 

57,5K^(.,)(.^) -:(zj ■ 

( Pj) (P2> — (P„) 

rhust be much less than 1. This is needed so that the steady state 
error to a ramp input will be as small as possible. The optimum 
case will be for one of z's = 0. 



3. 4 Compensation Schemes Attempted in the Feedback Path 
Around the Magamp. 

A compensator having the transfer function 



G(s) 



K(s^+2§ U) s + 
c z nz 

{s ■+ p^) (s + p^) 



O) ^ ) 
nz 



(3-7) 



was investigated. It was determined that if s = - 10 + j 10 were 
points through which it was desired to pass a 180 degree locus 
and it p^ = p^ = 230 were selected that using d pair of complex 
conjugate zeros ats = - l6 + j 10 would force a 180^ locus to pass 
through the desired points. 
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Since: 



57.5 K ‘ 

c 

Pi Pz 



«i 



57.5 K (356) 
c 



(230) 



« 1 



K << 2.585 

c 



Letting 1 does not satisfy specifications, but it is a 



(3-S) 



convenient place to start. Using Equation 3-7 and Equation 3-2 
yields: 



(s + 16 + j 10) (s + 16 - j 10) 
(s + 230) (s + 230) 



11.5 X 10^ s (s + 100) ^ (s^ + 10s + 100) 



(3-9) 



(s + 197) (s+135) (S+60+J43) (s+60-j43) (s +21+j28.5) (s+21-j28.5) 



An Esiac investigation performed on the root locus of Equation 
3-9 to find the roots of the characteristic equation with this com- 
pensator resulted in the root locas of Fig. 3-2. For K = 2 the 

c 

3 

system was unstable. By making K =6.17x10 the system 
could be made stable, but this violates Equation 3-8. This com- 
pensator was therefore adjudged unsatisfactory. 

A compensator having the transfer function 

(s + z) 



G(s) = c 



8^'+ lOs +100 ( s + P) 



(3-10) 



was investigated. 
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Substituting Equation 3-10 into Equation 3-2 .yields: 



K (s + z) 
c 



(s + 10s + 100) (s+p) 



- ,2 2 . 
ll,500(s) (s +100) (s +10S+100I 



(s+197) (s+135) (s+60+j43) (s+60-j43) {s+21+j2S.5) (s+21-j2S.5) 



(3-11) 



K ( s + z) 
c 



(3 + p) 



0 



11,500 s (s + 100)' 



(s+197) (s+135) (s+60+j43) (s+60-j4.3) (s+21+j28.5) (s+21-j28.5) 



= - 1 



By examining the pole - zero array of Equation 3-11 it was 
determined that z = 20, p = 183 would cause a 180° locus to pass 
through s = - 10 + j 10„ Substituting these values for z and p 



57. 3. K 



100 



( 20 ) 

( 18 ^ 



«1 



K << 16.1 
c 



(3-12) 



Again =1 does not satisfy this criterion adequately, but it 
make a convenient starting point. An Esiac investigation (Fig. 3-3) 
was performed on the root locus of Equation 3-11 and was found to 
be unstable for = 1.0. Stability could be achieved by setting 
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o 









K ^Z85, but this violates Equatiori 3-8. Therefore rfa? com- 
c i 



pensator was adjudged unsatisfactory. 



/ 



A compensator having the transfer function 

\ 2 



K 



G(s) = 



S -I- 2 



S + p 



s n+ 10s + 100 
was investigated. 

Substituting Equation 3-13 into Equation 3-2 yields. 

4 

,2 



(3-13) 



K ( s + z) 
c 



2 2 
(s +10s+i00) ( s+ 1 ) 



9 ") 

11,500 (s) (s + 10s + 100) (s t100)‘' 



( s+ 1 35)( s+ 1 9 7 ) ( s+ 60+ j 4 3) ( s+ 60 - j 43) ( s+ 2 1+ j 2S . 5) ( a+ 2 1 - j 28 5) 



= - 1 



K (s +z) 
c 



(3-14) 



(b+p) 



' h 



(11.500) (s + 100) s 



(s+135) (s+197) (s+60+j^3) (s+60-j43) (s+2i+J28,5) (s+21-j28.5) 



By examining the pole zero arr- ',y of Equation 3-14, it was determined 
that z =43, p =250 would cause a 180 locus to* pass through the points^,.,<l^ 






s = 10_+jl0. Applying the criterion of Equation 3-8 yields^ 






57.5 K 



100 



43 , 
250 



(3-15) 



K 58,7 

c 
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Assume = 1 satisfies this criterion. Performing an Esiac 
investigation (Fig. 3.4) on-the root locus of Equation 3-14 showed 
that the system was unstable for = 1 . Stability could be . • 
achi«v«d for ^ 432 1 but this viola,t«s Equation 3»15i Therefore 
this compensator was adjudged unacceptable. 

t 

1 ' 

,3. 5 Ross-Warren Method Used for Compensation 

From the root locus plot, it was desired to bring the unstable 

roots (s = 21.66 + 28. 6j) over into the left hand plane at a desirable 

location of s = - 10^ 10 j. This location should give a damping ratio, 

, within the specifications. This location should also give a 

good transient response. In order to obtain the original of 905, 

a 4 section compensator was tried. From the root locus plot, the 

following were obtained: 

0.= 46° 

4 

' ? = 56.35 

4 

Vk/G = .908 

cot 0, = .966 
4 

CSC 0 . = 1 . 39 
4 

^ = cot' ^ (cot 0^ - \/k/ G CSC 0^) 

= cot"^ (.966 - 1.26) = cot (-. 294 ) = 106° 
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■ TT + ®4 + S'- 



X 



= 180 + 46 - 56.35 - 106 = 63. 65 



and since ©^^0^, this compensator is allowable . Hence: 

z = h Vk/ G _ J = 18.908 sin 106 

^ sin 0, -J 



sin 63 . 65 



= 17,4 



P = 



h sin 



X 



4 ■ = 18 X . 956 



sin (0^- 0^) 



. 304 
56.2 



and the compensated system becomes 
- 1 



/s + 17.4 


4 


~11 .500 (s+100)^ s (s^ + 105 + 100) ' 1 


(s + 56.2 / 




(s+200) (s+100)^ s (s^ + 10s +100) +18; 1x10^°) 



This was tested for stability on the Esiac and the results are 
shown on Fig. 3.5. From thi it can be seen that the system was 
unstable (having two roots in the right hand plane). 

The next compensator investigated by this method was a 3 section 
filter. From root locus plot, 

0j = 6l . 3° cot 0^ = . ''47 

N^K/G = .88 . CSC 03 =1.14 

h = 18 
?= 56.35 
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+ + 







1 3 

X 3 = cot (cot 0^- \^KlQ 



CSC Oj) 



s cot"^ (.547-1.012^ 

= cot’^ (-.463) = 114.8 ^ 

= 7T+ ^ - X 3 - 180 + 61.3 - 56.35 - 114.8 

= 70.15° 

» 

and © 3 >' C5| ^ which shows that 3 sections is allowable. 



Hence z= h K/G 



sin )^ = 18 X .88 X .907 = 15.4 






sin 0, 



.940 



= h sin Xj 
sin (©2 -0^) 



= 18 X .907 



= 106 



.154 



The resulted in the following: 



-1 s+15.4 ^ ^ 11,500 (s-l-100)^ s (s^ + 10s + 100) 



s+106 



(s+ZOO) (s+100)^ s (s^+10s+100)+l8. 1 x 10^° 



This was tested for stability on the Esiac, and the results 
shown on Fig. 3.6 show that the system is still unstable (having 
two roots in the right hand plane s=26 + 36j). 



3.6 Compensator with Zero at the Origin 

A corhpensator having a transfer function 

r' t \ ^ s 

G (s) = c 

• s + p 



( 3 - 




>«H AND1»M 

lelsl t lal^ togM n 







was investigated. Substituting Equation 3-l6 into Equation 3-Z 

\ 

yields 



il.iOO i (8 + 100 )^ lOs + 100 ) K 8 

G 



(3-17) 



(s +197) (s+135) (S+60+J43) (s+60-j43) (s-21+j28.5) (s-21-j28.5) (s+p) 

t 

and it. follows that: 



57.5 K (0) 
■ e 






K oo 

c 

This satisfies the criterion for any finite value of K . The effects 

c 

of varying p were studied by selecting several values of p and plotting 
the root loci on the Esiac. It was determined that for values of /p/ — , 
4 , root loci 'from the two right half plane poles of Equation 3-17 

f 

came into ca,ncel the zeros at the origin; giving either roots in the 
right half plane or roots in the left half plane with real parts 1; 
neither of which is acceptable. 

The Esiac investigation further revealed that an optimum value 
of p was p^- 3.5. With this pole location, there could be produced 
two pairs of complex conjugate roots, the real part of each being 
= '^3.5. K required to obtain these root was K = 8. *96. The Esiac 
plot of this is Fig. 3.7. . 

However, the Esiac is not sensitive enough to determine the 

• ♦ 

actual value of K and root locations. Therefore a digital computer 
c 

program was utilized to determine the exact root location for various 
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values of K^. It was found that optimum results were obtained 
with K = 4.35. This gave roots as follows: 



s a . 2.46 + J 4.46 

s = - ’2 . 58 + j 5 . 925 

i 

s = - 101.904 
s = - 97.978 
8 = - 5.02 X 10^ 




These results satisfy all steady- state and transient specifications 
except settling time and the response to unit load torque. Checking , 
the reeiponse to a unit load torque, the block diagram of' the system 

is given in Fig. 3.8. The characterisitic determinant may be 

» 

written by inspection. (^^.3 ~ • 



A = 



1 

300iJ 



■ 50 1 + 

0 
0 



5xl0^s 



(s+200)(s+3.5) 



11,500 



0 . 0 



0 



8 + 200 



0 - 



10 



(s+lOO) 

0 



0 



9 43 

1 ■ . 0 
s+ . 246 



i2Ai 1 0 

8+10 

0 1 
^ S+.246 ^ 
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77 ^ 3 • 



For the unit-load torque input: 

ec s 



18 



r 1 
L ^ J 



8 ^0 



Obtaining each term separately: 



A ^ = ( 1 + 7.5s) 41 + 

s ^ o 



36,800 



( 3 - 18 ) 



f0c 



= 1 



s ^ 0 300s 



A 



ef 



= 36.9 + 2640s 



s^ 0 



Substituting these into Equation 3-18 



1 



ss 



©c = (36,9 + 2640s) 300s 

/, r V J, 36.800 

(1 + 71 . 5s) 41 + 



36.9 



(300)(36,800) 



8^0 



1 



300,000 



Radian 



And this is 



©c = .0115 minutes of arc, 
. ss 



3>7 Combining of Amplidyne Path with Tach Path 

Prior to attempting any compensation schemes an investigation 
was made to determine if there were any restrictions on the types 
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of compensators that could be used in these paths. Assuming 
an unspecified compensator in each path, the block diagt'aui of 
Fig. 3.9 is obtained. 



From th6 block diagram of Fig. 3.9, the characteristic 
determinant may be written down by inspection. 



(3-19) 



A= 



- 50 



315 X 10 



s(s+100)^(s^+10s+100) 



11,500 : 315 X 10 Kt Gc ,. 

s + 200 cl 100)^(s^+10s+100) 



11,500 



s + 200 



Now for a step input at 0 



ss 



step 



R 

^ac *^c0c 



s-0 



(3-20) 



Evaluating each term separately: 



A = 

8>0 



•50 1+57. 5G 



c 1 



■57.5 



. 315 



315K G - 
t c3 



A 905 

A =1 + 57. 5G, +18.1 K G.,+”^ 

s>0 t c3 




I 



41 










4 - 2 . . 



iCr 




= 2877 



A.. 



B>0 



°C0C 



.315 



8> 0 , S 
) 

Substituting into Equation 3-20 



E 



ss 



step 



1 - 



. 315 

2877 { s ) 



905 



1 +57. 5G , + 18.1 KG, + s 
cl t c3 



8 -> 0 



*^step 



1 - 



905 



57.5 G , 8 + 18. 1 K G- s + 905 
cl t c3 



8 -»> 0 



B.y inspection, it is seen that if it is desired that E be zero, 

s s 

step 



then G' , and G . can contain no terms of the form 1 
cl c3 — 

s 

^sl and the feedbacks be of opposite sign such that 

57.5 G., =- 18.11 K G , 
cl t c 3 



If however. 



■3.18 G , = K G 
cl . t c3 



Then the desired zero error will be attained. 
Now , 



1 



ss 



RAMP 



1 - 



A G 

ac C0C 



s ^ 0 



, (3-21) 



43 



I 



Substituting the already obtained values into this equation 



B 



ss 



RAMP 



. 901/ s 

1 + 57.5 G , + 18. 1 K G + 905/s 

cl t C3 



8>0 



_ 1 _ 

8 



1 + 57.5 G , + 18.1 KG 

cl t c3 

1 + 57.5 G - '+ 18.1 K G' + 905/ s 
cl t c 3 



_ 8 > 0 



1 + 57.5 G + 18. 1 K G , 
cj t ■ c3 

57.5 G , 8 + 18.1 KG ,s+905 
cl t c3 



8 ^0 



From this equation, it can be seen that G , and G , cannot contain 

cl c 3 

terms of the form 1 where N — 1 because this would cause the ' 

8 

numerator to go to oo . This then negates this possibility of a 

I 

solution using positive feedback in one path and negative in the other 
while using a term 1 as suggested from the equations for E 



ss 



step 



It is further seen from this equation that since the required 
E 



ss 



RAMP 



1 

905 



That the Lim G = Lim G , £ 0 

• 8>0 cl s > 0 c 3 



Considering now the possibility that the response to a unit load 

torque may place additional restrictions on the type of compensator 

adopted, use the block diagram of Fig. 3.8 with 4.35 s being 

8“+ 3.5 

replaced by a general compensator G^ ^ . 
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0c 



88 



u * 1 • t « 



A 



ef 



G 

£0c 



A 



8 > 0 



(3-22) 



Writing down A by inspection of Fig. 3.8 with 4.35 a ^ _ 

m :? = °ci 



A ' 
8-0 



0 



300s 



-50 1+57. 5G 



1 



0 G' 

t c3 

300 



■57.? 1 



100 1 38.3 0 



0 -1.04 1 



0 -36.9 1 

905 



= 1+57.5 G , + 18. 1 K G , +, 

cl t c3 s 



A , = 36.9 + 2064 G 

^e£ - cl 

8 >0 ' 



T0C 



1 

3uirr 



8-0 



Substityiting these equations into Equation 3-22, 












©c. 



1/300S [ 36.9 +2064 Gpl 



88 



U. 1. t. 



1 + 57.5 G ., + 18. 1 KG', + 

cl . t 63 8 



905 



8>0 



36.9 + 2064 G 



cL 



300s (1+57.5 G , + 18.1 K G- ) + 271,500 
• cl t c3 
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Now it can be seen that if the previously derived restrictions 






are observed; i.e. Lim G ' ^ 









tben ©c 



ss 



ult 



36.9 

271,500 



1 



7530 



c3 

t 

radian = . 468 minutes . 



Now that restrictions on the type of compensation have been 

r 1 

determined, it is useful to put the equations of ^ into a convenient 
form for working with G^^ and simultaneously. It is impossible 

to obtain a form of the characteristic equation in which G^j and 
G^^ appear simultaneously as cascade compensators. It is possible 
however to obtain an equation of a root locus which will be of the 
form 



n 



- 1 = 



K 7T (s + z.) G 
1 i=l ^ 



c3 



n 



■ ^2 jil' 



m 



n 



I I (s'+ p ) 
k=l 



or qonversely: 



1 - 



KjTT 



n 



1 + (s + z.) G , 

4j.l ^ j' c3 



n 



Now clearing tjie denominator: 



4 = 



1 + 



18. lxiO^®(8+lO0)^9*+lOs+lOO) G , + 36.2xlO*K 80 . 

' cl t c3 

s + 4 1 0^+5.41x10^8^2.54x10^5^2.5x10^8^+2x10^3+18. 1x10^® 



Letting this 6th order equation =, 



18. 1x10^°(b+100)^{s^+10s+100) K N , 

A= 1 + 






D 



cl 



36.2xlO®KsN , 
t c3 



D 



A 



D 



c3 



where 



G , = K , N 
cl cl cl 



D 



cl 



G = N 

c3 c3 



D 



c3 



Now iA may be put into one of the two previously mentioned forms. 
The first form is : 



A. 



1 + 



36. 2x10 K sN - 
t c3 

°A 



s(s+100)(s +10s+100)N , K , 

cl cl 



1 + 



cl 
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Clearing the denominator 



36.2x10® K 8 N .D , 
4 . 1 1 



^ 1 ^ 1 0 8+ 1 00) K^jN^j J 



(3-23) 



The bracketed term in the denominator may now be. used to form 
the root locus equation 



s(8+100)^(s +I0.+ 100) K ,N , 

A = 1+ £i_£i. 



(3-24) 



Thus it is seen that the effect of the compensator is to supply . 

✓ 

zeros in Equation 3-23 for each pole of G and also to provide 

c'l 

the poles of Equation 3-23 in the forrh of the roots of Equation 3-24. 

The alternate form of the equd’tion, by the principle of similarity 
is: 



. 18.1x10^® (s+100)^(s^+10s+100) D - N ,K , 

^ A D + 36.2x10® K s N -) D , 

D c3 t c3 cl 



Of which the denominator bracketed term may be taken to form 
the root locus equation 



A = 
8 



36.2x10® K s N , 
t c3 

“2SI D~" 

D c3 
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Assviming a compensator G 



4.35 G 



cl 8+3.5 ’ 



This function for is the compensator which almost 

m^t specifications previously. It is therefore taken as a 

starting point to find a Q , which will actually meet specifications , 

c-5 

Substituting this equation for G^^ into Equation 3-Z3: 




1 + 



36.2 X 10^ K s (s+3.5) G , 
t c3 

Roots obtained in Section VII 



(e) 



It is seen that the pole zero array of Fig. 3. 10 will result. 

Now G ^ must be a compensator containing a term (s) in the,^'^’’^ 

.numerator. This will result in a double zero at the origin of 

Fig. 3.10. Obviously, a root locus will always appear in the region. 

between one pair of the complex poles at s =-2.5 ahd the origin. 

Since this will always give a set of roots with their real parts -^-2.5 , 

this compensator will not produce satisfactory results in terms of 

tran^eht response settling time specifications. 

It can be seen that to achieve satisfactory compensation with 

the assurqed G , , it will probably be necessary. for G , to cancel 
cl c-^ , 

all four complex poles at s = - 2.5. If the compensator 

be this complicated, it would be just as well, to use a complicated 

» 

compensator for G , and achieve compensation by use of the mag 

cl 

amp path; alone , 
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It may be seen further from the sketch of Fig. 3.11, the 
root locus of 

*A = I’i 1 Ij 5008(8+ i 00) (b" + 10s+100) _ 

A ~ G^l (3.26) 

D 

where it is required that G contain a term (s) in the numerator 

cl 

that th^re will always be at least two roots of Equation 3-26 in the 

I • 

area between s=20 and s = - 5; since the root locus from the right 
half plane poles will always go to either the zeros at the origin or 
at 8 = - 5+ j 8.6. These roots would then become poles of the 
equation 




36.2 X 10 K s D , 
t cl 



n 

TT 

i=l 



c3 



(s+pi) 



and since as mentioned previously contain a term s 

in the numerator there would in all probability be root locus between 
the poles under discussion and the double zero at the origin; unless 
of course, produced zeros to cancel those poles. Since it is 

necessary to retain the zeros at the origin, the obvious alternative 
is to cancel the zeros at s = - 5 + j8.6 with an active compensator 
fqr . If this be done, all specifications can be met using the 
mag amp path only and no additional advantages accrue from also 
compensating with the tach path. 
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The denominator of each of these root locus equations is 
another root locus. The denominator root locus is a function of 
only one compensator thus one compensator may be used to deter- 
mine the poles which the other compensator will work with. In 
many cases, this may be a useful tool. Obtaining these forms 
for the two compensation paths now under consideration: 
Expanding the determinant of Equation 3-19i 



^=1 



18.1 X 10 



10 



s(s+100)^(s+200)(s^+10s+100 






s + 200 



cl 



36.2x10 KG _ 
t c3 



(s+100)^(s+200)(s^+10s+100) 



11,500 



= 1 + 



s+200 cl i 



1 + 



18. lx 10 



10 



s(s+T00)^(s +200)(s^.+ 10s+100) 



. 8 . 



36.2x10 K G - 
t c3 



(s-HOO)^ {s+200)(s^-H0s+100) 



1 + 



18.1 X 10 



10 



s(s+100)^s+200)(s^+10s+100) 
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When using an active compensator in the rhag amp path to 
cancel the zeros at s = - 5 -I- j8.6, ^s a first consideration this 
compensator must be cascaded with a compensator of the form 

Q ’ * 

g ^ p • An active compensator may be used to produce poles ' 



at s = - 5 +' j8.6. ‘There will be two real zeros which come as 
an integral part of this compensator,' which may be placed as 

a 

desired. The farther out on the real axis these zeros are p^ced, 
the lower the root locus gain •will become to place roots in the . 

I. • 

desired locations. This presents no real problem since a voltage 
divider network may be made part of the compensator and only a 



f 



part of the mag amp output feedback as compensation. 

M I ■ . 

Therefore, for simplicity, let the two real zeros cancel 
the poles at s = - 197 and s = - 135. This ghres the equation* 



I 



ll,500s(s+100)^(s^+10s+100)K ,s (s+135)(s+197) 

;r — 2 

(s+p) (s flOs+100) 



+ 1 



n.500 ‘^ci . 

(s-21+j28.5)(s-211j28.5)(s+60+j43)(s+60-j43)(s+p) 



+ 1 
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It was determined from Esiac investigation that p = - 3 would 
produce a root locus through the desired ^ range. (See Fig. 3. 12) 
The root locus gain to produce roots with real part =7.5 was 102. 

« ' I 

The compensator then become s 



,, (=^>35) (.+ 197). 

= — 2 — — 

^ (s + 10s + 100) (s+3) 

, 4 

Substituting this into Equation 3-26 and setting = 0, the 
characterisitic equation polynomial becomes: 

0 = s^ + 181.68s'^ +22,035.85s^ + 1 , 095 , 320 . 975^+8 , 506, 777 . 62s+ 
25,070,121 

Solving this equation on the digital computer gave the following 
roots: 

s = - 4.47 4 2.88j 
s = - 53.28+ 103j 
s = - 66.18' 

Investigation of these roots showed that the required specifica- 
tions were not met. Several values of gain were used until satisfactory 
roots were obtained to meet specifications. The value of K used 
was 39 . 0 giving the following roots: 

V 

s = - 42.6+ 60j 
s = - 1 1 . 2 + 20j 
s = - 16.0 
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These roots were obtained from the^ root locus plot shown in 
Fig. 3, .13. Since the compensator satisfies the steady state, 
erlteirla and ihdse roots satisfy the transient specifications, this 
compensator was selected. From Fig, 3.13, it can be seen that 
the dominant (s = l6.0) root is far enough away from the origin 
to give a settling time within the re^quired specifications. A 
damping coefficient of S' = .485 resulted from the above roots. 
Since the inclusion of the tach path offered no advantages in 
simplification of the compensator required, it was not investigated 
further. Analyzing the transient response gives; 



^c _ 10.1 X lO^^s (s-l-8) (s^+ 103 -I- 100) 

. (s+16) (s+11.2 + 20j) (s +42.6 + 60j) 

From this, it can be seen that the settling time is 4 x _ = 

'1 1 « ^ 

i 357 seconds which is well within the required specifications . 

I 

This detailed illustration is intended to show some of the uses 
of the determinant-block diagram relation. The solution obtained 
is not a complete one , nor does the problem itself have any bearing 
on the self-adaptive problem other than as an illustration of manipu 
lative methods. 
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4. 


DEFINITION OF SYMBOLS USED IN FOLLOWING SECTIONS 


U 


A 


x component of velocity (scalar) 




V, 


4 


y component of velocity (scalar) 




w 


^ m 


ii emnpamm 




A 


A 


2 

Moment of inertia about x axis, slug-ft 


(whole mass) 


B 


A 


2 

moment of inertia about y axis, slug-ft 


(whole muss) 



4 
A 

4 

4 
4 
4 
A 
A 

K, = 



C 

D 

E 

F 

P 

Q 

R 

K, 



a 

b 



moment of inertia about z axis, slug-it (whole mass) 

2 

product of inertia, j yz dm, slug-ft 

product of inertia, ^ xz dm, aiug-fi^ 

^ 2 

product of inertia, j xy dm, slug- ft 

angular component of rotation about x-axis (scalar) radian/ sec 
angular component of rotation about y-axis (scalar) radian/ sec 
angular component of rotation about z-axis (scalar) radian/ sec 
relative angular velocity of rotor term on x-axis (scalar) 
relative angular velocity of rotor term on y-axis (scalar) 



K- = relative angular velpcity of rotor term on z-axis (scalar) 



4 
4 
A 

c = 

A 

a = 
c 

m ^ 
m — 

A 

e = 



moment of inertia of' rotor along x-axis ^ 

moment, of inertia of rotor along y-axis 
moment of inertia of rotor along z-axis 

2 

acceleration of airframe at its mass center ft/ sec 
mass of control surface 

eccentricity of contrbl surface between hinge and the control 
mass center 



0 = 



angular displacement about the reference (x^) axis 
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t 

I 

r 

I, 

ii^ 

i 



e ^ 
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4ir rp 
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4ii ''r 

r 

/^Uo^SrCr C 

:r. h 

2.r rj, 

/’UoSrCr C, ' 

4xr 

the elevator position angle (radians) 

h 

the sileron position angle (radians) 

II i' 

the rudder position angle (radians) 
force in the direction of the x axis (scalar) 
force in the direction of the y axis (scalar) 

force in the direction of the z axis (scalai^ 

) 

moment about the x axis (scalar) 
moment about the y axis (scalar) 
moment about the z axis (scalar) 
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5. DERIVATION OF THE STANDARD BLOCK DIAGRAM FOR THE 

AIRFRAME. 

5. 1 Introduction 

The types o£ motion that result from a disturbance in some 
equilibrium flight condition and the transient motion of the aircraft 
in response to control movements are analyzed in the study of the 
dynamic characteristics of an airplane. Over the course of years the 
effects of dynamic stability in aircraft have been known to designers. 

Tj'he mathematical computations were discussed by early aeronautics 
pioneers which included Lanchester , Bryant and Glauert. In spite of the 
fact that dynamic stability effects were known, most people paid little 
concern to this field. As design groups solved other problems of air- 
craft design, the dynamic stability effects were automatically solved as 
a secondary result. However, since World War II, aircraft development 
has pushed into the field of very high performance aircraft. In many 
cases the dynamic stability problem doesn't solve itself as happened 
previously. Furthermore, the fire coritrol problem has entered the, scene 
which has generated a control' problem of no small magnitude. 

This work will consider the aircraft in six degrees of freedom and 
explore paths of compensation and soliition. It is the prime intention to 
explpre the work of Chu (1) in multi loop servo systems and determine the 
value of it as applied to the aircraft stability. ' 

The first section will set down a general definition of terms and 

. ( 

basic equations of motion as related to aircraft. It will further set,^'^ 

down the basic assumptions uf^qn which 'the equations are formed. 
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The equ.aticoG motion for the airframe which are completely 
derived in Appendix I are: 



= (U + QW - RV) m s X - mg Bin e 


(5-1) 


. 1 

7 {V + RU - PW) m = Y + mg cos 0 sin V 

« 

F = {,W + PV - QU) m = Z + mg cos e cos Y 

1, = AP + aKj + (C - B) QR + D (Q^ - R^) - E (PQ .+ R) 


(5-2),, 

( 5 I 3 ) 

1 


+ FPR + cK^Q - b K^R 

4 

M = BQ + bK^ + (A - C) PQ + E (P^ - -R^) - FQR + DPQ 


(5-4) 


+ aKjR - cK^R 

1 

N = CR + cK^ + (B - A) PQ + F (Q^ - P^) + E (QR - P) 


(5-5) 

, » 


- DPQ + bK^P - aKjQ 


(5-6) 


H+F=It{rmla =P (RP - Q) 
e e e *' e e cz ex 


(5.-7) 


H+F=Ih =mla -(R+ PQ) P - (RQ - P) rz 
r r T9 r r cy rx 


(5-8) 


2H + F =I X +2P (RQ + P ) 
a a a 0 ay 

The angular and linear velocity of a firee body in space are: 
P = 0 - ''V sin 0 • ■ 


(5-9) 

(5-10) 



(5-10) 



(5-11) 



» 

Q = © fos © + V sin 0 cos © 

R = V cos 0 cos 0-0 sin 0 (5-12) 

-r~ - U cos © cos V + V (sin 0 sin © cos V - cos 0 sin V ) 
at ' 

1 + W (cos 0 sin © cos V + sin 0 sin V ) (5-13) 

* U cos © sin + V (sin 0 sin © sin + cos 0 cos V ) 
at 

+ W (cos 0 sin © sin - sin 0 cos (5-14) 

-j2 — - - U sin © + V sin 0 cos © + W cos 0 cos © (5-15) 

at 

Consideration of the aerodynamic forces must be resolved into the 
effects they produce on various parts of the body in some flight condi- 
tions and also the resultant moment effects they produce. 

The first or necessary condition is to establish a reference condi- 
tion and proceed from there. A reference condition is some flight 
condition where equilibrium exists in all equations and the summation 
of forces is zero. This is a steady state condition. 

, ^ I 

The set of axis upon which this rigid body is oriented is commonly 
called the Eulerian Axis system. General aerodynamics usage modifies 
the z axis so the positive direction points downward. 
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The first rotation moves about the Z axis in a direction IP and is known 
as yaw. The next is the movement about the Y axis which is term © 
and known as pitch. The third motion about an axis is roll about the X 
axis and is maneuvered in terms of 0, Caution must be used in deter- 
mining these positions because the order in using this system is not 

\ 

commutative . 

The derivation for the yaw rate R, the roll rate P and the pitch 

4 

rate Q are shown in detail in Appendix Two. They are: 



P = 0 - y sin 0 (5-16) 

\ 

Q = 0 cos 0 + lf> sin 0 cos © ' • ' (5-17) 

R = ^ cos 0 cos 0-0 sih 0 ' (5-18) 

and where small angles are involved 

P'^ 0 !, '1 (5-19) 

© (,5-20) 

R^ V I' . (5-21) 



Similarly» Etkin shows: * 

0 = Q cos 0 - R sin 0 
0 = P + Q sin 0 tan © + R cos 0 tan © 
= (Q sin 0 + R cos 0) slec © 



(5-22) 

(^- 23 ) 

(5-24) 



The general equations of motion have to be applied to the aircraft motion. 

I 

The effects of axis movement' upon the'^weight force of the aircraft will^*”^ 
be shown as well as the reference state condition. It must be r'^membered 
that the gravity force is always oriented toward the center of the<earth. 
The' reference body axis may or may ilot be oriented in this direction. 
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NORMAL AIRFRAME AXIS ORIENTATION 
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From the discussion of the Eulerian angle concept, it is deter- 
mined that a reasonable reference condition would be 
F = O = X - mg sin 
> V - W eos sin 0 

O O O t 

F s 0 = Z .cos 0 cos 0 
z o o o 

t 

This allows for a steady state condition where the aircraft is not 
aecejerating but neither is it necessarily in straight and level flight. 

4 

From the previous illustration of Eulerian angles coupled with the deri- 
vations in Appendix II, the gra'vity forces can be determined as shown 



be low . 


,1 






r * * 


cos © cos Ip 


cos © sinlp 


-sin © 




my 

X 


cos V sin 0 

/ 


cos (pcos 0 


cos © sin 0 


xmg= 


my 

y 


- sin cos 0 


+ sin ^ sin © sin 0 








cos sin © cos 0 


sin 4^ sin © ctis 0 


cos © cos 0 




my 

z 


- sin ^ .sin 0 


-cos ^ sin 0 









(5-25) 

(S.26) 

(5-27)-' 



jC: r.cafi sion and derivation of this gjeneral cosine matrix are 
found in several sources (2,3). The overall result of the equations 
of motion is shown below. ^ 

« 

’ * i-, ’ 

. .... _ ^ I 

. • . !' 
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. m (U + QW, - RV) = .X - (mg sin B^) (cos 0 cos V ) + mg (cos 0^ sin 0^) 

(c/bs © sin V ) - mg (cos 0 cos '0 ) (cos © sin 0 ) •. (5-28) 

o o 

I • 



> m ( V + RU - PW) a Y • mg (sin 9^) (cos 4^ sin O sin 0 • sin V oos 0) 

+mg (cos © sin 0 ) (cos cos 0 + sin ^ sin 0 sin 0) 

■ 0 .0 

•+mg (cos © cos 0 ) (cos © sin 0 ) (5-29) 

o o 



0 

m (W + PV - QU) = Z - mg ( sin ©^) (cos 4^ sin © cos 0 - sinV sin 0) 

+ mg (cos © sin 0 ) (sin V sin © cos 0 - cos 4^ sin 0) 
o o 

+ mg (cos © cos 0 ) (cos © cos 0 ) , (5-30) 

o o 



L = AP + aK^ + QR (C-B) -E (PQ + R) + cK^Q - bK^R 

• ' .• 2 i’ ’ 

M = BQ + bK^ + PQ.(A-C) + E (P - R ) + aK^R - cK^P 
N = CR + cKj + PQ (B-A) + E (QR = P) + bK^P - aK^Q 



(5-31) 

I 

(5-32) 

(5-33) 



. I 

The above equations are complete fexcept for the external forces 
X, Y, and Z and the external moments’ L, 'M, and N. These forces and 
moments are composed of aerodynamic and propulsion effects and the 

moments due to control surfaces. These equations are very general in 

I 

nature and contain the following assumptions for simplification from the 

most general and complex case. These assumptions are: 

\ 

1. The earth. is assumed to be fixed in space and the Earth's 

, 

atmosphere is assumed to be fixed with respect to the earth* 

2. The airframe is assumed to be 'a rigid body. '' 
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3. The mass of the body is assximed constant during the time 

’ of the problem. 

4. For an aircraft or missile, the XZ plane is assumed to have 
mirrer symmetry. 

The above equations then cover the six degrees of freedom found in an 
aircraft or missile. Coupled with the other auxiliary equations covering 
roll, pitch and yaw effects along with control motions, solution is pro- 
bably possible. To date, this has been accomplished either by a machine 
using specific coefficients for various terms or several simplifications. 
Discussion of this will be conducted in succeeding sections. 



i 



I 
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5.2 



Conaideration of External Forces and Moments on a Body in 
Motion 

In the previous section the equations of motion were discussed. 

The aerodynamic thrust and control effects were simply grouped in X, 

Y. Z, Li, M and N. It is the aim of this chapter to discuss these in 

detail with respect to the effect on the dynamic control stability problem. 

Consider the term X when the aircraft is in a disturbed condition. 

The reference condition is X and thus X' = X + ^X for some disturbed 

o o 

condition. The question is theh to determine the ^ X components. 

The force and moment variables can be expressed in a Taylor's 



Series . 

F = Fo + ( -f5 -) 



^0( 



2 






0 ^ 2 ! 



(5-34) 



i : .1 

For practical reasons the effects of second and higher order terms are 

I • 1 1 1 ■ I ' 

omitted from further considerations as will be discussed later. Con- 



i' ■ 'I 

sidering the motion of the aircraft itself the X forces would be 



(5-35) 



X = xo ) u t ) u ) W + w . (— ) q t (-^) , 

••I I 1 

+ (“T ) V + (— r-) V + ( ) p + ( J-) p + (— ) r + (-^-T- ) r 



Similarly, 



(5-36) 









^ W 






^q' 






For convenience* all terms will be expressed in letter subscript form>*"" 

,1 

It can readily be seen then that considering X, Y, Z, L, M and N and 



i I 
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the effects Qf U, V, W, P, Q, and R plus the acceleration terms U, 

V, W , P, Q, and R there are 7Z derivative terms required. This does 
not include the qontrol effects which can add another. 36 derivatives. 
Needless to say» considerable effort has been expended in justifying 
simplifying assumptions. In addition to this, obtaining reliable co- ' 
efficients for various derivative terms is not always possible. 

I 

It is not the intention of this work to cover the computational 
methods for these various derivatives. The scope is too large. A 
Search of all authoritative books j^3, 4, 5, 6, Jon the subject yields a 
table of terms which can be reasonably determined. This appears as 
Table 1. These considerations are mentioned to bring out the awareness 
of the different causes and effects of external forces. There are con- 
siderable gaps in the table. Determination of these to complete the 
table would be a topic of a complete study in itself. Much work remains 
to be done in this field. 

I Furthermore, some of the methods of computing the knovirn 
stability derivatives’are questibnable . Some computations from theory 
prove quite accurate. Others are very poor and empirical methods Or 
model testing are used. Even these fail at times and "educated guesses" 
from previous airplanes are used as a basis of determining stability 

derivatives for a new design. Tail effects are the most notable in this 

♦ 

last category. ' 
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TABLE 5-1 KNOWN AERODYNAMIC STABILITY DERIVATIVES 
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In all aerodynamics problems, considerable effort is made to 
non-dimensionalize the equations. The advantages of doing this are 
of considerable merit. For the most part, it is rather simply done. 

1 % been i^und in ^ein iheery and aetuai praetiee that 

F=l/2/^V^SC_ (5-37) 

r 



Where S is a reference area, usually the wing area, V is the absolute 
reference velocity and is some non.dimensional coefficient. C^, 
varies with Mach Number and Reynolds Number depending on the force 
considered, (lift, drag, etc.) ^ is the density of the fluid through 
which the body is passing. Thus a term such as Zw has the following 
relationship: 



Zw' = 



^ w 



= 1/2 /^U 2^ (-^^) =1/2 /'U ' SC = 1/2/^SU C 
f o S ^ V/ /oZ f o 

u w 



(5-38) 



w 



Here again V and S are the velqcity and wing area terms respectiveljy . 
The same concept applies to the moment terms except that a moment 
arrnternq must also appear in^ the equa^^ion. A generalized example is: 



M 



= Pxr SIC 
/ 0 m 



(5-39) 



In this example M is some moment and C. , is a needed dimensionless 

M 

coefficient to satisfy the equation. ^ and S are the same as 

mentioned above while ^is the^ moment arm. 

Likewise in a stability moment the derivative can be of the following 



form? 

I. 
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(5-40) 



= SU^l(-^-^)h= U' SjC 

h o 3 h - o m. 

Another external force to be considered is that of thrust. 

This bsGomss important factor in most flight conditions. , The. 

location of the thrust axis is seldom coincident with any reference 

axis in the body. 

n , * 




FIGURE 5-2 

' 1 



I 



AXIS SYSTEM FOR LOCATING THRUST IN AN AIRFRAME 
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Assuming the thrust line lies in the xz plane, then the thrust 



forces can be resolved in this generalized manner. 



= T cos e .J, 

= T sin © .J, 

= T 

T z 

i 

For a steady state condition the equations become 
X = To cos ©„ 



Z = To sin OT 
°T 



M = T 
Orp oz 



<5-41) 

(S-42) 

(5-43) 

(5-44) 

(5-45) 

(5-46) 



while in a disturbed condition where the body axis remains fixed in 
tbe mass, 



X.„ = T, cos ©^ 
T 1 T 



(5-47). 

, ■^i'5-48) 

, (5-49) 
(5-50) 

Assuming the air density remains constant so the, engine thrust will 



Zrj, = Tj sin 

M = T 
T Iz 



and Tj = To + A T 



not be affected 

4 

3 T 



^ ^ <^RPM 



(5-51) 



or 



T=Tu+Tw ^RPM 

u ^RPM 



77 



These are the considerations for the external forces on a body 

4 . * 

in motion. While the dimensionless coefficient has been used exten- 
sively by the aeronautical engineering field* it will not be used further 
in this work. The primary mission of this work is to investigate ^e 
control response and paths for possible compensation'. , 

. I 

The generalized equations of motion then appear: . ,, 



• • 

. m(U + QW - RV) + frig sin (cos © cos.M* ) - mg (cos ©^ sin 0 ^) 

(cos © sinV) + mg (cos © cos 0 ) (sin ©) = X + X* u + X* w + 

® o o o u w 



(5-52) 



X' - w + X'^q + X'^r\ + i + X\v^'T^ cos ©^ + T^ cos ©^u + 

1 ^ rpm 



w q*’. 'I' ■‘■v' ‘o T ' ‘u T 

C? rpm 



T cos © , ^ 



(5-53) 

e 

m(y ■+ RU - PW) ‘+ mg sin ©^ (cos V sin © sin 0 - sin ^ cos 0) - ,mg 

(cos © 8in 0 ) (cos V cos 0 + sin Vsii' © Bin 0) - mg (cos © cos 0 ) 

O O . 0,0 

» 

(cos © sin ©) = Y + Y' V + Y' .V + Y' P + Y' r + Y'* X + Y' ? + Y'/lT 

o . V V p^ r, S 

\ (51-54) 

m(W + PV - QU) + mg sin ©^ (cos 4^ sin © cos 0 + sin V sin 0) 

- mg (cos © sin © ) ( sin ^ sin © cos 0 cos sin.©) - mg . " 

^ o o 

(cos © cos 0 ) (cos © cos 0) = Z + Z' u + Z' w + Z'* w + Z' q + 

• o o • o u w w q’ 

e 

Z'T*n. + Z'^f\, - T sin ©_ - T sin ©_u j T w__. , sin ©_ 

H o T u T ^ ^RPM T^RPM 



(5-55) 



AP + aK, + QR (C-B) - E (PQ + R) +cK,Q - bK.,R = L + L' v + 
1 • I 3 2 o V 

+ L'pP + L'^r + C. 
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(5-56) 



BQ + bK + PQ ^A-C) + E(P^ - R^) + a K R = cK P = M + M' u 
it X o VX 

+ M' w+M'.w+M' q+M'.q + M'^H. + + T z + T u 

w ,w q q ^ ^ o u 

\ 

' (5-57) 

CR + cK, PQ (B-A) + E (QR-P) + bK^P = aK O = N + N' v + 

3 2 1 o V 

N'pP + N'^r + N'^ J + N'^jr + 

Thus, it has been shown that the equations of motion have the above 
external forces which will affect the rigid mass in motion. The assump 
tions made are: 

1 . The earth is assumed to be fixed in space and the Earth's 

atmosphere is assumed to be fixed with respect to the 
I Earth. 

2. The mass is a rigid b|ody. 

3. The mass of the body is assumed constant during the time ' 

of the problem. i| i' 

i 4. For an aircraft or missile » the xz plane is assumed to have i 
mirror symmetry. 

These assumptions were made in the first section. In acldition the 
following assumptions and simplifications are made in considering the 
aei;odynamic forces : i 

I ' ■ 
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1. The flow is assumed to be quasi- steady . 

2. The effects of auxiliary aerodynamic equipment such as 
flaps I speed brakes, slats and spoilers are neglected. ,, 

3. The thrust can be affected only by change of speed and/or 

revolutions per minute. No Change is considered for change 
in fuel mixture, manifold pressure, side slip or propeller 
fin effect. * 

These assumptions were made for several reasons. In the first and 
third cases, absence of adequate information on many of the coefficients 
makes ignoring or setting them equal to zero mandatory. In all probability 
these coefficients are very small. Assximption two covers special cases. 
While these items may have considerable effect when used, they are 
eliminated from the general ciise at this time because an aircraft pro- 
bably uses them less than two 'percent of its flying time. For the special 
cases, these effects may be quickly added. 

The work ir> these first two sections is developed primarily for 
background purposes. There is nothing new or startling in this. For 
further examination of this subject, the reader is referred to the works 
in Ref. 7 and 8 which closely parallel the development in this work. 

The aerodynamic force and moment terms carry a primed notation 
for ease in publication. In lat^r sections these terms will be div^ided'by 
mass or moment of inertia. These terms after division will not carry 
the primeC ) notation. There are far more of the latter terms used. Thus, 

it is for ease of publication only and does not affect theory in any way. 
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5.3 The Airframe Equations in a Standard Form 



This section will take the general equations developed in the 
previous chapters and simplify them. The simplification will extend 
to an aircraft in straight and level flight and subject to small pertur- 
bations . This is the standard method of dealing with the problem as 
found in present day texts « However, these equations will be used 
in block diagram form and in determinantal arrays as developed in 
Chu*s work. (1) The results will be compared with present day 
classical methods to determine validity of this determinantal concept. 

The first simplification is the concept of the lateral system being 
independent of the longitudinal system. Therefore, the six equations 
of motion are divided into the two groups. The longitudina] group 
includes the directional motion in the X and Z directions and rotary 
motion about Y axis. The lateral group consists of directional motion 
along the Y axis and rotary motion about the X and Z axis. The axis 
system referred to here is the body axis system of the airplane. 

Secondly, the aircraft is initially in a steady state condition. 

This is the reference state as mentioned previously. When the dis- 
turbance terms are all zero Equation 5-52 reduces to: 

m (Uo + Qo Wo - RoVo) + mg sin ©o= Xo 4- To cos ©^ (5-58) 

However, it was stated that this situation was restricted to longi- 
tudinal or lateral systems. Thus for motion along the X body axis, 
a yaw motion Ro and a side motion Vo would be zero. Further it was 
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stated that the aircraft was in a steady state reference condition. 

9 

Therefore U would be zero. Thus Equation 5-58 reduces to 
m (Qo Wo) + mg sin = X'o + To cos © T (5-59) 

Similarly, it may be observed that Equation 5-53 for the steady 
state reduces to 

Y'o = 0 (5-60) 

Likewise Equation 5-54 reduces to 

m (Qo Vo) - Z* + mg cos © 4- To sin 0 T = 0 (5-6l) 

o o 

The moment equations all reduce to zero provided the rotor effects 
are neglected or considered zero. 

Lo = Mo = No = 0 (5-62) 

Recalling that only small disturbances were to be considered 
and these disturbances were to be in only the symmetrical or 
longitudinal system for the longitudinal group the derivative due 
to assymmetric variables of V, P, R, ^ and ^ are zero. In a 
similar manner, when the asymmetric system effects are considered 
the symmetric disturbances are ignored. This leaves the variables 
V, W, Q, and )^out of the asymmetric equations. Finally, because 
of the small disturbances, the change angle © will be small so 
the approximation cos © ^ 1 and sin 0*= © is valid. 

The result of the previous discussion is subtracting Equations 
5-59 and 5-52, eliminating the asymmetric changes and supplying 
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t’he small disturbance considerations. The result is: 



m(u +) +mg e cos © =X “U + X »W + X JW +X «Q 

on w w q (5-63) 

+ X^T| + X^«)^ + Tu' cos G^U + T^« ®T 

For the z direction taking 5-61 from 5-54 and using the same 



assumptions: 

m (w-) + mg ©^ sin 0 =Z«U+Z"W+ZJW+Z*Q+ (5-64) 

o u w w q ' ' 

# 

Zy\'^ + Z • ' 0 - TixU sin - T o ^ 

L ^ ^ ^ T ORPM ORPM sin ©T 

For the moment in the y axis, the following applies: 

BQ =M^U + MwW + MwW + MqQ + MqQ + + (5 65) 



T ’U - + T <5'rPMz 

u zl 5RPM 

To complete the case for the symmetrical group, the action of the 
elevators must be included. From appendix I the elevator hinge 
moment effects are reduced to : 






H n 

e>^ 



i- HewW + 



HewW + 



I 



(5-66) 



HeuU + HeqQ + ^F e 



In most cases where moment of inertia terms are handled it 
is usually expedient to use an axis system which is the principal 
inertia axis. This means that the product of inertia terms are zero. 
For aerodynamic motiorx problems, this applies for some cases. (8) 
This case is not one of them. What is gained in simplification of 
the rotary motion equations by use of principal axis is lost in 
attempting to handle the force equations with the wind velocity 
considerations in sine and cosine terms. By orienting the X axis 
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(5-67) 



into the reference wind velocity the result is: 

2 2 2 2 
V = U + V + w 
c 

where U = U + u. 

o 

Therefore, in the steady state reference condition V = U and 

c o 

V = W =0. With the assumption that u, v and w were very 
o o 

small, their products and squares can be neglected. 
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V \/ U ^ + 2U u 
c ^ o o 

Z 

and since U is much greater than 2 U u, then V = U = U and these 
o o c o 

terms are interchangeable in this section. 

As previously stated, the reference situation is a steady state 

motion and includes first order effects only. The orientation of the 

axis U is a finite value while V = W =0. Furthermore, the author 
o o o 

is hard pressed to visualize an aircraft or missile tumbling along with 

a motion although P exists in spinning missiles and bullets*. Furthei 

X =X = Z • = Xyy = = =0as an assumption . 

w q w 'I ^ (^RPM ^ 

These terms are usually very small. 

The ’’standard" determinant form as stated in Chapter 3 of Chu^ is: 



1 + PF^ ^^11 


PF^ 


PF 

21 


PF 

3 


PF 

31 


-PF. 


1 + PF^ 


PF 

22 


PF 3 


PF 

32 


-^^13 


- 




1 + PF 


PF 

33 



where the PF'»^ indicates a performance function. The system may have 
n nodes. This example contains only three nodes. The "standard" 
block diagram for this form is: 



^•'In sections 2 and 3 the symbol G was used to designate the transfer func- 
tion of a component, and the component itself was normally a distinct piece 
of hardware. The performance function PF is also a transfer function, 
but as used with regard to the airframe it designates a relationship be- 
tween two variables which are not visualized as the input and output of a 
piece of hardware. The use of PF notation from here on is therefore a 
deliberate choice fo emphasize this distinction. 
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FIGURE 5. 3 STANDARD BLOCK DIAGRAM DRAWN FROM 3x3 
DETERMINANT. 
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The nodes are labeled in Roman Numerals. The basic rules 



of nodal signals allow several inputs, but only one output. This 
output from the node goes directly to a performance function with a 
single subscript number. Signals may be picked off from the node 
output and fed forward to another node. An example of this is con- 
tained in the signal going through PF^^. In this case the signal is 
also multiplied by the function PF and then fed to node III. Feedback 
takes place in a similar manner. The signal output from the first system 
PFj is picked off and fed through back to node I. Similarly the 

other nodes function in the same manner. 

For feedback blocks, the first subscript number denotes the out- 
put signal source for the feed back and the second subscript number 
the input node to which the feedback signal is going. The notable 
points about the determinant are: 

1. The main diagonal always has terms of one plus the 

direct path performance function times the self-loop 
feedback performance function. 

(1 + PF^PFjj) 

2. All terms above and to the right of the main diagonal 
are feedback terms* 

3* All terms below and to the left of the main diagonal of 
the ‘^standard" determinant form are of a feed forward 
nature. All of these feed forward terms are negative. 

The term immediately below the main diagonal term is the 
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direct path performance function. Thus, where the second 
column main diagonal spot has 1 + ^^22' direct path 

term will appear in the location immediately below as 

-PF^. 

These rules are strictly mechanical. For a rigorous derivation of this 
system the reader is referred directly to Chu’ s work. 

In Chapter 9 of Etkin (4), there exists a two degree of freedom 
situation where © and W vary and the speed in the forward direction is 
constant. This is the short period approximation. The equations using 
the LaPlace transforms and nondimensional coefficients are: 



(2^ s - G 


) OC - {2^ + G^^) s 0 . Gz^-Y^ = 0 


(5-68) 


-G OC 

mot 


2 - ^ 

+ {i^s = G )© - Gm T) =0 

p mqs 1 ^ 


(5-69) 


-°hoc^ 


-°h, 


(5-70) 



s ^ e 



This author also grouped his acceleration terms such that G, 

included C, and C, , . Similar situations exist for G, , G 

hyn hoc moC 

and G . Note that r = ^ and P = 0 for this consideration. From 

m^ 

these equations a block diagram is evolved. 
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FIGURE 5.4 BLOCK DIAGRAM OF TWO DEGREES OF FREEDOM IN 
THE LONGITUDINAL SYSTEM SHOWN BY ETKIN^ 

Further, Etkin derives: 






s - G^^) + G 



zn 



moC 



(5-71) 



’'l (2^ S - J , 



Jcjt 



Ac(. 



G ( i A s - C ) + G ( 2^ + G ) 
Zyi K mq m Zq 






(5-72) 



n. yi 



<2>< s - Gza:> - °mq ' - °Zq> * 



Now it is intended to use Chu*s work involving the use o£ '^standard” 
block diagrams and ’’standard” determinant form to prove that the 
’’standard” system applies for the airframe as well as a multiloop servo 
system » 
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The Z, M and equations of motion with the applied assump- 



tions and © = O are: 
o 




• 

w-Uq-Z W-ZQ - Z^^ - zVi^= 0 
o w q *L ^ 


(5-73) 


• • ' 

&-M W-M.W-MQ - Mys"^- M)o 71=0 

w w q H ^ 


(5-74) 


vi- lOn'*^ lOyi^- 10 W - 10 .W - 10 q=AFe 
^ 1 ^ w w q 


(5-75) 



Recalling that 60= Uo tan OC and for small angles 
^CCand 

The equations in LaPlace notation then become: 

(s - sio 10 (sio^ + 10^) (slOq) e = A Fe (5-76) 

(s^ - sM e - (sMy^ + My^)^ - (sMC< + M OC ) Uo Ak =0 (5-77) 

(s -Z^) U^AoC - (sU^ + sZ^) e -n (sZ>j^+ Z )^) '>^ = 0 (5-78) 



These equations in standard determinant form now become: 



n 


A OC 


e 


, (slOiq + 10>x ) 


(sl0<x + lOOC) 


slOq 


^ ■ 2 


2 


2 


s 


s 


s 


(sZ/) + Zy] ) 


, zo^ 


s(Uo + Zq) 


— V 

s 


s 


s 


( sM y[ + M y 7 ) 


(sMci+ MO() 


^ sMq 


2 ‘ 


2 


^ ■ 2 


s 


s 


1 

1 
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Ordinarily, it is not necessary to divide by s or s However, 
it IS done here to demonstrate more closely the standard formo 

To take the differential equations and arrange them by rows and 
columns to see the actual physical system is to no avail. The arrange- 
ment does not exist in the above form. It will be shown how to achieve 
this form later. 

Consider the ^standard form” array as set forth above. The 
block diagram for this in Chu's "standard form” is: 




FIGURE 5.5 RESULT OF BLOCK DIAGRAMING EQUATIONS 5-76, 
5-77, and 5-78. 
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The above block diagram then is a three node system but the outputs 
from the first and second main path blocks are nothing that is readily 
recognized in an airframe. Because of the fact that the main block 
following the third node is 1 , the output here is ©. 

A few simple block diagram manipulations clears the situation 
and shows true signal outputs that are physically realized in the air- 
frame. Block diagram manipulations allow the changing of a pick off 
point provided the transfer function is corrected. Further, two blocks 
can be combined into one equivalent transfer function. 




which 

changes 

to 




Figure 5.6 EQUIVALENT BLOCK DIAGRAM TERMS 



Adding a dummy node this now changes to: 




Figure 5v7 FURTHER BLOCK DIAGRAM COMBINATIONS 
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The result of using these two maneuvers rearranges the block diagram 



to the following: 




Figure 5.8 REARRANGED BLOCK DIAGRAM OF FIGURE 5.5, 
USING MANIPULATIONS OF FIGURES 5.6 and 5.7. 

Maneuvering the above diagram location, the block diagram takes the 
form of Fig. 5.9, which coincides with Fig. 5.4. 
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Figure 5.9 TWO DEGREES OF FREEDOM IN AN AIRFRAME 
LONGITUDINAL SYSTEM . 

(rearranged diagram of Fig, 5.8 to compare with Fig. 5.4) 
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It IS noted that this block diagram approaches the form set 



down in Etkin except for nondimensionalizing the terms. The dashed 



area in the diagram is what Etkin labels G_ and Gotw. 

The determinant from the above dashed portion of the block diagram 



is: 



1 


0 


0 


0 


0 


(sZri + Zy^) 


j ZOC 
s 


0 


-s(Uo + Zq) 
s 


0 


0 


- 1 


1 


0 


0 


-{sM')j + M^) 
s ^ 


0 


- (sMO( +M0() 
2 


1 - sMq 

2 

s 


0 


0 


0 


0 1 


-1 


1 



0 ^ oC 

Evaluating for the transfer functions and — — , the results 

, . Gjki!\iy 

will be in the form 

This form is the type expressed by Chu (1) and Dialer (12), The 
term, which is the whole determinant, reduces to: 



A = 



1 - Zoi 
s 

-1 

0 



1 

-(sM«^ + MOC) 
2 



• s(U^ + Zq) 



0 

1 - sM 



which is = (1 - ZqC ) ( 1 - sM ) - s (U + Z ) (sM^ + MC^) (5-79) 

, q Q q 

e c 2 Q c 2 



To consider the term the cofactor Aiy of this determinant is 

obtained by crossing through the first row and the second column. 



-A 



12 = 



-s (U + Z ) 

o q 



s 

0 



s 

0 



( sM -t- M y| ) -( sM OC -1- M ) 1 - sM^ 



12 = +'l 
3 



+ (sZ + Z r| ) (s - sM ) + s (U^ + Z ) (sMy^+Mr|) 



The main feed forward transfer function between the input Signal from 
node 2 to the output is +1 in this case. This is also true of node 4 
input to output. Hence = + 1* 

The transfer function 



'AOCY[ 



Ad: 






= 



12 



A 



(5-81) 



G^ocr| = 



(sZy|^ +Zy|)(s - sM ) + s { U o + Z ) { sM + M y| ) 

(s - Za) (s^ - sM ) - s(U + z ) (sMOC + MCX ) 

q o q 



It can be noted that this equation 5-80 is the same equation as 
5-72 given by Etkin save the non-dimensionalizing terms. 
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To evaluate the transfer function © 1 , the same method is used 



'^46 = ■ ‘ 



A 



14 = 



-{sZ-r)+ Z 






1 - ZoC 



■(sMyj_4_M_^) 0 



-(sM<^ + M<X) 0 



A 



14 = 



- (sZ + Z y| ) 1 - ZOC 






-1 



.(sMCX + MCK ) 



A 



(5-82) 



14 =_1 (sM>^+My^)(s - ZQ() + ( sMOC + M0() (sZV|^+ Zy^) 



Therefore: 



*^40 A 4 

A 



^ _ (sM>^ + Mn ) (s - Zcc ) + (sMO, + MqQ (sZ>^ + Zy^) (5-83) 

(s - ZCX) (s^ - sM ) - s (U + Z ) (sm6(+ MoC) 

q o q 
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Here again, excepting the non-dimensionalzing, the result is the 
same. Equation 5-83 is equal to equation 5-71 as given by Etkin, 

The system of ’’standard” forms then works for the airframe 
differential equations after some manipulation. It can be noted that 
there is a method of making up this standard form without the necessity 
of the manipulation which will yield physically realizable output signals. 
When the block diagram of Fig, 5.8 is converted into a determinant 
for solution of the transfer function, there appears a series of alterna- 
ting rows where the equation -1 + 1 = 0 appears. It also appears that 
these simple equations occur in a sequence where the plus one term is 
placed in the 2n ! row and column, n is the number of differential 
equations in the problem. The minus one term is placed in the 2n ! 
rows and (2n !-l) columns. All feedforward terms (those below and 
to the left of the main diagonal) are moved one column to the right 
and remain in the same row. This doubles the size of the determinant. 
However, there are several z^ro terms introduced and the actual work 
involved in evaluation of the determinant is not markedly increased. 

The next item is to look at the lateral motion problem. It is 
slightly more complex because there are five differential equations 
of motion instead of four. 

Equation 5-35 with simplifications becomes: 

(v + Uor) - 0 mg cos 0^ - YvV - Yp - Y - Y^^ = 0 (5-84) 
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(5-85) 



( A + r) - (J) g coa© 
' o 

U 

o 




Y^r - Y 

■u“ U 

o o 



= 0 



The moment equation 5-55 becomes 



P R 
A 



LpP- V - 






= 0 



(5-86) 



and equation 5-57 becomes 



R + ^ 
C 



^ - N^/5 -N^P - N^R - ^ 



(5-87) 



The aileron hinge motion equations from Appendix I are 



UpP- ll^r- -n§& = A 



Fa 



(5-88) 



And also from Appendix I, the rudder hinge motion is 



^ -12 12pP - 12^r - 12^^ - 12^5 = AFr (5-89) 

From these terms a determinantal array can be made. It should 



be noted that this time a set of equations of - 1 + 1 = 0 will be inserted 



in the system. This permits direct drawing of the block diagram and 



omitting the need of manipulation. Table 5-2 then presents this array. 



From this Fig. 5.10 is drawn directly. 
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DETERMINANT ARRAY OF LATERAL SYSTEM EQUATIONS FOR SMALL DISTURBANCES 




lOI 





) The resulting block diagram very clearly parallels that of the 
Si:eral system shown by Etkin in Chapter 9 (4). The "standard" form 






k achieve this block diagram differs 

1 There is, then, the question of why does this differ from the 
j'iginal "standard" form set forth by Chu. 

i 

Chapter Two of Chu (1) discusses at some length the idea of 
irangment of varying performance functions into a "standard" block 

'f 

diagram. The individual pieces of hardware are at hand with known 

.dividual performance functions. The orderly arrangement of these 

prmit the construction of a ’’standard” determinant for the whole system. 

From this, a multiloop servo system can be handled in generalized 

1 

results. In the case of the airframe, the determinant formed by the 

1 

differential equations was known and the block diagram had to be 

V 

‘irmed. Admittedly a block diagram (Fig. 5.5) was made up but the 
:itput signals are not the physically realizable type that are generally 
lown to pilots and engineers. The same performance function results 
om the system. Both are mathematically correct. 

The difference lies in the fact that Chu is taking individual pieces 
' hardware where the functions are known. The blocks are made up 
id then the ’’standard” determinant form is made. In the case of the 
Lrplane, the equations of motion are known and a determinant is 
nmediately formed^ Each column contains the performance functions 
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multiplied by the direct function signal output term. Chu’s standard 



form dictates that all terms above the main diagonal are multiplied by 
the node output or input to the direct function. Steps must be taken 
then to change this in the differential equations. The addition of the 
alternating -1 + 1 = 0 equations solves this very nicely and in no way 
alters the value of the equation^ 

A series of evaluations of determinant forms is given in Appendix 
III, The results show that the addition of these terms in no way changes 
the final value of a determinant or a cofactor value. 

Referring to Fig. 5.8 there are other block forms that may change 
the values within the blocks, but will not alter the final outcome. Con- 
sider the first node in Fig. 5.8 as an example. The four forms expressed 
in Fig . 5.11 are equivalent . The original form is shown in Fig . 5.11a. 

The self-loop feedback path is combined with the direct path in Fig. 

2 

5.11b. The s term is eliminated and the performance functions are 
then in the form shown in Fig. 5.11c. A further change can be made by 
dividing the nodal input terms by the denominator of the direct path 
feed-forward function. The result is then in the form shown in Fig. 

5. lid. The performance functions in each block differ somewhat, but 
the final result remains the same. These manipulations allow the 
user more freedom in the forms of the performance functions. It 
allows changes in form to fit the needs of the particular problem. 
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Figure 5.11 (a) Figure 5.11 (b) 




Figure 5.11 (c) Figure 5,11 (d) 



FIGURE 5.11 PERMISSIBLE BLOCK VALUE VARIATIONS AT A 
NODE. 
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The pick-off point for the feed forward terms does differ in 
requirements. It is seen in Figs. 5.11a and (d) that the pick-off can 
be either off the output signal of node I or II. The final result will 
not change . However, this is not the case in Figs . 5. 11b and (c). The 
pick- off point must be off the output signal of node II. The final 
result will differ if the output of node I was tapped in the latter case. 

For converting a set of linear differential equations into the 
standard block diagram form, the following rules can be applied. 

1. Between each differential equation, an equation -1 + 1 = 0 
shall be inserted. An additional equation -1 + 1 = 0 shall 
be added at the bottom of the set so the number of additional 
equations shall equal the number of differential equations 
originally in the problem. 

2. The arrangement of the values of these added equations in 
the determinant shall be such that the positive one value 

is always in the main diagonal block. The minus one term 
shall always be placed in the block immediately to the left 
of the main diagonal term. 

3. All other blocks in the rows of these added equations shall 
be of zero value. 

4. All finite performance functions of the original differential 
equations which are located below the main diagonal shall 
be moved one place to the right. The vacated blocks shall 
be replaced with a zero value. This last rule is optional 
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to the user. It does save the possibility of error which 
could appear because of cases shown in Fig. 5.11b and (d). 

With the simplified rules of making a standard block diagram 
from the standard determinant, a valuable tool is developed for many 
purposes. It can provide a student with a means to draw and visualize 
a complex feedback system such as an airframe from a group of differen- 
tial equations. For the research groups, it makes analog computer 
programming a simplified matter with respect to amplifier usage. In 
this latter case, pitfalls do exist. The airframe is no simple matter. 

The reader is referred to Howe (8) for an excellent coverage of the air- 
frame problem for analog computers. Most important, though, is the 
fact that this standard form of Chu's is valid for any multi- loop linear 
feedback system whether it is electrical, mechanical or a combination 
of the two. Thisthenputs a regular form to the multi-loop control 
system from which design or synthesis can proceed in an orderly manner. 
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6, Expansion of the Standard Form for Greater Disturbances 

The system developed by Chu was for linear functions. It is 
the purpose of this section to explore the equations of motion and 
determine whether the equations of motion can be used for larger dis- 
turbances than considered in Section V. If so, then to what degree 
can this system be used. 

For all practical purposes, the small angle approximation holds 
quite well up to one -tenth radian. Therefore, sin © and 

cos 0^1. With this simplification in mind, Equation 5-52 is then 

rewritten. 

* 

m ( U + QW - RV) + mg sin © - ip mg (cos © sin 0 ) + 

o' o o 

• 

© mg (cos © cos 0 ) = X + X' U + X 'W + X'.W + X' Q + X'Y,fl 
o oouww q 

t X^>i + X'^V + cos cos + T j JRPM cos 

(6-1) 

The thrust forces are combined into the forward velocity motion 



u and a steady state reference condition of 
m (Q W - R V ) + mg sin © = X + T cos ©^ 



o o o o 
is assumed 



Subtracting 6-2 from 6-1 the following equation results 

m (u + Q + W + qw - R - V r - vr) 
ow oq ov o 

- y^mg (cos © sin 0 ) + ©mg (cos © cos 0 ) - X* 



- X* - X’. W - X»Q - X 



w 



w 






x^''! - x'v = 0 



( 6 - 2 ) 



(6-3) 
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The product term qw is dropped as a small value product situation. 



The same applies to the vr term 

Converting to La Place notation and dividing by the mass the 

following equation will result 

+ e U sW e - Lp U £ - Ip g cos Q sin 0 

00 O ^oor o o 



su 



+ © g cos © cos 0 -XU- (sX • -f X ) ^ - sXq © 

o o u w w 



(sx-;i^+x,^)r|^ - = 0 



(6-4) 



The y and z direction force equations will form up in the same manner. 

sU .+^u+U - W 0 - tpg sin e 

00 o o o ° o 



g (cos © sin 0 ) - g 0 cos © cos 0 + XaA 

o o „ o p 

o ' 



(6-5) 



sY 0 - sY^ (f. (Y^ - (sY^ + Yg')'? = 0 



sU AC< + P U /? + V 0-Q -V 0+gsin©© 

O O O ^ os OU OS o 

+ g cos © sin 0 0 - g cos © cos 0 - Z U 

o o o o 

- (sZw + Zw) .^0^-sZ © - (sZ'Y|^+ Z ^ = 0 



( 6 - 6 ) 
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The moment equations can be formed in much the same manner. From 



Section Five, Equation 5-55 would convert to this equation: 

) (sQ \V + sR e ) . ^ (sP e - sQ ip 

A o' o A o o A ^ 

(6-7) 

- (sL^ + )^- 0 - (sL^ 7^)-(sL^ +L^)^+L^^ 

In a similar manner Equations 5-56 and 5-57 can be brought to a similar 
maximum linear condition. They are shown in the following equations: 

s^© + ) (sP 0 + sQ 0) + (sP 0 - sRo ip ) 

Boo ^ o 



-Mu 

u 



( sMw + 



M ) 
w 



U ^0C-s(sM*+M ) © 

n n n / 



(sM 






+ 



M 7^ ) T[ = 0 



(6-8) 



m + (sP © + sQ 0 ) + (sQ ip + sR © - s^0 ) 

^'C o o C o 

- -sNp 0 - Ip - N|^ ■ +N^)? (6-9) 

Also the control hinge moment Equations 5-66, 5-88, and 5-89 
are recalled in their same form. 
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These nine equations represent the linear equations of motion 
of the airframe and include all cross-coupling effects These equations 
still retain the requirement of a linear approximation The equations 
with the author's modification of adding the equations -1+1 = 0 make 
up Table 6-1 . This is an 18 x 18 determinant From this table the 
block diagram of Fig. 6.1 can be drawn. 
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TABLE 6" I Determinant Array of Airframe with Linear Approximations for Six Degrees of Freedom. 
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This block diagram then represents the complete airframe 
response maintaining linear approximations. Even then this system 
has limitations. As an example, consider some airplane flying with 
a steady state wing angle of attack of 15 degrees, (assume a ’’con- 
vention” wing and the reference axis is the zero lift line) Any pilot 
will quickly testify that an increase of 5 degrees angle of attack at 
that point will probably result in a rather wildly oscillating maneuver 
popularly known as a spin This condition is generally outside of the 
range of the aerodynamic forces and moments usually computed. 

Some aerodynamic coefficients such as lift (C ) are decidedly non- 

JL< 

linear in this region. Therefore, the restrictions imposed on this 
derivation of terms must include a restriction that the airframe is 
not operating at the very edge of its flight envelope. 

The system has been diagrammed directly from Table 6-1 . The 
terms in Column One are all connected at node I of Fig. 6- 1 . Likewise, 
the terms in Column Two of Table 6-1 are all connected to node II. 

Extra nodes appear in the beginning and end of the block diagram. 

These are inserted as provisions for control functions. These control 
nodes are lettered alphabetically. The block diagram contains the 
required loops for aileron and elevator control. The rudder control 
loop can be easily inserted when desired. 

Examination of Fig. 6.1 reveals many interesting points here- 
tofore known but sometimes only by intuitive reasoning. The concept 
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of separate systems for the symmetric axis and asymmetric axis 
systems is very well shown. The level flight steady state conditions 
used in Section Three stated that 

Po = Qo = Ro = Vo = Wo = 0 (6-10) 

Examining the paths between nodes VIII and IX reveal that all 
of the performance functions are zero under these specified conditions 
of steady state flight. There is no direct, feedforward or feedback 
paths between nodes VIII and IX. Thus the restricted case of steady 
state level flight conditions with small perturbation is clearly divided 
into the two systems. The longitudinal system includes the forward 
velocity, vertical velocity, pitching motion and elevator position. 

The block diagram for the system can be further restricted to two 
degrees of freedom by holding the X velocity change, u, at zero. 

This allows change in vertical motion and pitch. This would then be 
Figs. 5-8 and 5-9 of Section 5. 

The so-called lateral airframe response system comprises the 
nodes from IX to XVIII. The rudder and aileron motion, sideslip, 
yaw and roll comprise the motion for this system. With the steady 
state level flight conditions previously mentioned, the block diagram 
is identical to Fig. 5-10 of Section 5. 

An interesting note of these separate systems is the cross 
coupling terms that begin to arise with larger disturbances. An 
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important one is the pitch due to sideslip. The effects are labeled 

in the diagram. Several aircraft are subject to this. Usually, 

this is a tucking or nose down motion that appears as sideslip angle 

increases. In the simplified systems this effect was neglected. 

It is of negligible value for very small disturbances, but can become 

noticeable by the time the sideslip angle ( ^ ) reaches five degrees. 

The effect exists and is shown by this expanded system. 

Excellent methods of analysis are available for these standard 

forms brought forth in this work. Methods of analyzing outputs in 

response to one or more inputs are available. Work by Thaler (10) 

contains a very useable system. 

Where N = the output signal at the n th node 

^ = characteristic determinant 

A = the cofactor of the determinant with "a” 
n the node location receiving the input 

and ’’n" the node location of the output 
signal. 

A 

Gj = the performance function on the input 
^ signal prior to entering the input node. 

I = the input signal 

G = the direct path performance function 
nN 

between the output node and the output 
signal. 

Thus a single input signal would yield an output signal 

A G, I 
a 1 
n a 

L 



Nt = G T.T 
I nN 



( 6 - 11 ) 
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For two or more input signals additional subscripts would be used 
and the functions are additive. The rules of superposition of linear 
outputs justify this method. 



N = G 



nN 

A 



A G, I+Zl, G, II+A G- III+--- 
a 1 b 2., c 3 

n a n b n c 



(6-12) 



This equation would give the resulting ouput of one or more inputs 
at various nodes in the system. 

The degrees of freedom can be restricted also. In the block 
diagram, (Fig. 6-2) consider a cut made in the path at point "s". 




FIGURE 6-2 ILLUSTRATION OF RESTRICTING FREEDOM OF 
MOVEMENT IN A SYSTEM BY CUTTING A NODE OUTPUT IN 
THE BLOCK DIAGRAM. 
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This opening in the path now eliminates the signal to the 
direct path, the feedback and the feedforward path. The result 
of this in the determinant is the insertion of zeros in all the 
performance functions in the column. It leaves only the constant, 
one, in the main diagonal block. The reader will recall that 
the main diagonal has the term 1 + PF^PF^^. This is shown 
below as opening the output path from a third node. 
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The output path from node three is opened at point s and the deter- 
minant IS changed in this manner □ 

A special case occurs for the diagram illustrated in Figo 6-2* 

The second node has only one input signal. Because the path was cut 
at point ”s”, the output from the second node will also be zero. Thus, 
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all feedforward values from this node will be zero. The direct 
and feedback paths of the second node are non-existant to being with 
and remain so. The determinant form would appear as now shown. 

X X 0 0 X X 



X X 0 0 X 



X 



X X 1 0 X X 



X X 0 1 X X 



X X 0 0 X X 



X X 0 0 X X 



Simple manipulation now reduces the six by six determinant to a 
four by four determinant. Therefore, when it is desired to consider 
a problem with one or more degrees of freedom held at zero, the 
simple method of opening the output signal path from the node and 
placing a zero for all the performance functions in the column leaving 
only the constant value one in the main diagonal point accomplishes 
this . 

Referring to Section 5, the reader can observe the general 
equations of motion for mass particles and examine one such as 
Equation 5-55 which is repeated below, 

AP + a^k^ + QR (C-B) - E (PQ + R) 4- cK^Q - bK^R = 



+ L;v + L'.^V + L'^p + L'_r + L'^5 + ^ 



O V 



(6-13) 



+ - 0 
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Note the third term QR is a product of (Q + q) (R + r) . The pro- 

o o 

duct when multiplied isQR + Qr+Rq+qr. From the equations 

o o o o 

of Section 5: 



When the small angle approximation is no longer valid, the handling 
of this product term becomes rather involved. Furthermore, the aero- 
dynamic coefficients of the L' , L* , L’ etc terms lose their linear 

V p r 

approximations over large changes in flight conditions. Most aero- 
dynamic coefficients are of a decidedly non-linear nature over the 
full span of an aircraft's maneuvering capability. The coefficients 
are definitely non- linear with Mach number when the velocity range 
covers the subsonic, transonic and supersonic spectrum. The reader 
can observe the general equations of motion in Section 5. Considering 
the complexities of incorporating equations 6-14 and 6-15 for large 
movements and the non linearities of the aerodynamic coefficients, 
it is not long before one realizes that the non linearities are numerous 
and cover all the equations. The prospect of block diagramming for 
the purpose of placing all non linearities in one block to be handled 
by a describing function scheme is hopeless. 

The possibility of doing the problem with an analog computer 
and using no approximations is fair. Undoubtedly the task, if attempted. 



Q = © cos 0 + 




(6-14) 



J I • 

R = ^ cos 0 cos © - © sin 0 



m 



(6-15) 
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will be of considerable magnitude. 

Thus, it is seen that this system of standard determinant arrays 
and standard block diagrams can be used to handle the airf rameproblem 
as long as a linear approximation is maintained. The system is 
flexible and can carry varied conditions of multiple inputs. It can 
also be noted that for the six degrees of freedom problem over a large 
variation, it is of limited value except for diagramming a computer 
problem . 
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7. APPLICATION OF THE DETERMINANTAL APPROACH TO 
A SIMPLIFIED AIRFRAME CONTROL PROBLEMS. 

The system to be considered is a surface to air homing missile 

using proportional navigation, A sketch of the missile shape is 



guidance computer provides command signals for the yaw and 
pitch autopilots, which have the same form because of the symmetry 
of the missile. The missile is to be considered a rigid body, and 
gravity forces are considered to be trimmed out in steady flight. 

The missile is also to be considered in steady state operation 
such that the speed (Mach) is constant and the pitch angle of attack 
IS constant at OC » Assuming that cos ^^/^=1.0, sin 

assuming that the mass and moments of inertia are constant, the 
linearized equations of motion of the missile can be derived and 



shown in Fig. 7.1, which also defines the coordinate system. A 



are: 



Yaw force equation 




e 




(7-1) 



o 



o 



Yaw moment equation 






= 0 



(7-2) 



Roll moment equation 




(7-3) 



Where the coefficients are as defined in Section 4. 
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Fig. 7-1 sketch: of missile shape 
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Since pitch is constant, the equations for the roll and yaw 
control systems are needed. The measuring devices used are 
accelerometers which measure acceleration in the y and z 
direction or, within a scale factor 

and + Y 




Other measuring devices are rate gyros which measure 0, and ©. 
The relationships existing in the controlled missile are defined 
by the functional block diagram of Fig. 7.2. 

If the sensors are considered perfect G, =G., =G =G =G^= 1.0, 

12345 



and the control equations may be written; 

sji =(-) [ 



0 + 0 + K, 0 

4 5 6 






+K^ i\> +K^4^ (7-4) 



(7-5) 



The coefficients of these equations may be arranged in a deter- 
minantal array as shown on Table 7-l« 

A block diagram may be constructed directly from this ^ , 
then manipulated using the methods indicated in Section 5, or some 
preliminary algebraic manipulations may be used to define transfer 
functions. Using the latter scheme, each row is divided by the 
entry in the main diagonal of that row, providing a new determinantal 
array {Equation 7-7) shown as Table 7-Z and 7-2a. 

The standard block diagram may be drawn from Equation 7-7, 
and inspection of the diagram indicates that a rearrangement of 
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FIGURE 7-2 




Where i K , K_, K., K_ , K/ = gain constants G, , G^ , G_ , G., 
c 5 ^ o o 1234 



G , = transfer function of the form 

^ 2 

W 

G = ^ 

^ s^+2 J W s+W ^ 

n n 
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TABLE 7- 
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Equation (7-6) 



1 


o 


O 


o 


o 


1 

o 


O 


o 


o 


1 

o 


r-H 




> 


L/* 

>- 

1 

d 


O 


M/' 

2 

\ 




o 


t 


1 

d 


o 


» 


CL 

-J 

1 

cL 


o 


- 


7 


1 

1 


o 


o 


o 


o 


o 


o 


o 


f*H 


3^ 

'~'+ 


<i 

< 


O 




o 


o 


o 


o 


o 

'O 

c 

*io 

\ 


I 

Ci 


o 




r-H 

1 


o 


o 




o 


o 


o 


o 


o 


r-H 


1 


Q. 

1 

< 


o 


o 


o 




1 


1 

ci 


o 


z 

1 




o 


1*^ 


f-H 

1 


o 


o 


o 


o 




o 


o 


o 


i-H 


0 

^t5 

CO 

o 

-1- 


1 

<1 

Vw 


o 


o 

1 


o 


o 


o 




L 


1 


o 


- 


r-H 

1 


o 


o 


o 


o 


o 


o 




o 


r-H 


r/> 

z 

1 


1 


o 


Xo 

> 

1 


i 


o 


A/1 

_1 

1 


o. 
, -J 
1 

< 


o 


o 


o 


Co 




i-H 


o 


o 


o 


o 


o 


o 


o 


o 



19 .^ 



TABLE 7-2 Equation (7-7) 



Where the performance functions of Section 5 are defined as: 







(kj s + K^) 

(s - ) 

*^1 V 

(s-KjY. ) 
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(s-Lp) 




-cosc^ ^ 
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? 3 J (s-^> 




^ (s-Lp) 




jS 



(s-Lp) 



TABLE 7 -2a 
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nodes and blocks can produce a separation of the roll and yaw 



subsystems. If this is done the block diagram takes the standard 
form shown in Fig. 7.3. Note that the yaw and roll subsystems 
are clearly separated, and that specific blocks are obviously 
representative of cross-coupling effects. A load disturbance 
^ ^ is indicated, and a command input is also shown as 
(Y + Y ) c. The input block PFa is defined 



as 






-K 






= PFa 



)c 



5-K Y, 






It can be noted for a simple missile, various cross coupling 
terms are zero at .zero angle of attack in the asymmetric system 






N 



(s-K^ ) 



1 






(s-Y^ ) (s-Lp) 



= 0 (7-8) 



and since CC = 0, then sin o( 



o . . 1 II lii *- =0 

(s-^) 

Examining Fig. 7-3, and noting the zero blocks resulting from 
Equation 7-8, the standard determinant can be formed as shown 
in Equation 7-9- (Table 7-3) 
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Roll to yaw couplings 



Yaw Subsystem 

-K3S+K2 



■ N. 
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Roll subsystem 



TABLE 7-3 



Equation 7-9 
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Therefore, the selection of input nodes in the order of 
^ ^ and ^ has allowed the formation of a yaw 

and roll subsystem. The array clearly shows the two sub- 
systems and the two attendant cross coupling areas 

Specifications for dynamic response of simple missiles is 
often most rigid at the zero angle of attack condition. It can be 
noted from Equation 7-9 that compensation can be made separately 
in the tvo subsystems. After these specifications for zero angle 
of attack are met, the cross coupling consideration may be 
brought into the problem to check design requirement at other 
angles of attack. 

It is evident that this resultant block diagram is in a form 
which conveys a maximum amount of information to the designer. 
With proper manipulations the design of a well controlled aircraft 
or missile becomes a much simpler problem . 

No attempt is made here to apply the results just obtained. 

The purpose of this section was to indicate how the block diagram, 
determinant technique can be applied to the simplification of 
complex airframe control problems . It is felt that the preceding 
considerations have accomplished this. 
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8. THE GENERATION OF COMPLEX ZEROS 



8 o 1 Introduction 

Any closed loop system which can be represented by a block dia- 
gram has a characteristic equation which may always be manipulated 
into the form 



n 

K jr (s + z ) 

— G = . 1 (8-1) 

TT (s + p ) 

j=l ^ 



Where , 0 ^ K - oo < m < , and the z. and p may be real and/ or 

1 j 

complex, K is assumed to represent a variable gain element of the 

system. represents the transfer function of a compensator which may 

be in any single feed forward or feedback path, or may be cascaded in the 

main transmission channel. G is of the form 

c 



G 

c 



n 

K TT (s + z.) 

^ i= 1 



m IT 

j = 1 



(s + 






(8-2) 



The definitions of the terms in Equation 8-2 are the same as Equation 
8- 1 . The total gain of Equation 8 . 1 is then 



K_ =KK (8-3) 

RL c 

The selection of a transfer function for G , besides being based upon 

c 

physical realizability is guided by the desire to force the roots of the 
characteristic equation to lie in such locations on the s-plane that the 
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transient response of the system to be compensated will meet specifica- 
tions without affecting steady state accuracy. Steady state accuracy 
specifications are normally met before compensation by setting system 
gain to obtain an acceptable steady state error coefficient o However, 
it is conceivable that compensation could be required to meet steady state 
specifications . 

It is common practice to select two complex conjugate root loca- 
tions and then to choose G so that the two selected locations will be roots 

c 

of the characteristic equation. Generally, it is hoped that the selected 
roots will prove to be dominant, but a certain amount of experience and 
foresight is required in order to ensure that this end will be attained. 

Even with experience, some trial and error is frequently involved. In 
any case, only calculation of the transient response will determine whether 
or not a particular solution is acceptable. The concept of dominance 
states that if the residues at certain roots of the characteristic equation 
are large compared to the residues at all other roots, and if the real 
parts of those certain roots are several times smaller than the real parts 
of any other roots, then those certain roots are dominant ^ A transient 
response calculated using only the dominant roots will be sufficiently 
accurate for engineering purposes . All roots must be included however 
when calculating the phase angle of the transient response. Note that 
this concept may significantly reduce the labor involved in calculating 
the transient response of higher order systems. 
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The common use of passive networks to effect compensation is 
with containing only negative real poles and zeros, and usually with 
K^4^1. It IS possible to generate complex poles or zeros with passive 
networks, but it is difficult to use these networks practically due to large 
attenuation in the compensator and extreme sensitivity of the complex 
pole or zero location to variations of components within the network. In 
some cases an additional amplifier is required to correct for the attenua- 
tion of a passive compensator, so the compensator could have been de- 
signed as an active compensator from the start. 

The terms ’’active network” and ’’active compensator” are inter- 
changeable and refer to networks which contain passive components, a 
summer, and one or two amplifiers. 

Frequently, complex poles and/or complex zeros are present in a 
system before compensation is attempted. Having networks available to 
use as compensators which will produce complex poles or zeros will in 
some instance make the compensation problem easier. The active net- 
works investigated in this report are flexible in that the complex poles 
or zeros resulting from them may be relocated on the s-plane by varying 
either an amplifier gain or the magnitude of one or more of the passive 
elements in the network. 

The passive networks, which are essential to active compensators, 
may be either RC or RL networks. However, RL networks pose some 
problems due to difficulties in determining and allowing for the resis- 
tance of the inductive elements. 
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Active compensation is not envisioned as replacing passive com- 
pensation, but rather as supplementing it by providing solutions to 
compensation problems not easily attainable with passive compensators 
alone. In some instances it may prove desirable to combine active and 
passive compensators to achieve desired results- 

Two major disadvantages of active compensators are cost and com- 
plexity. Both of these factors arise from the fact that amplifiers and 
summers are required as parts of the compensator. 

The passive networks, their transfer functions, and s-plane pole- 
zero arrays which will be used as basic components in the active networks 
are depicted in Fig. 8.1. They are held to simple configurations pri- 
marily to provide attainable goals for this investigation. In addition, 
all compensation problems attempted could be solved using compensators 
which required only these simple networks. The transfer function of 
Fig. 8. 1 are defined on Table 8-1 for algebraic manipulation. 
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NETWORK 



S- PLANE 



TRANSFER FUNCTION 



(a) 



(b) 
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L R^+R^ L 





s 

s + R/L 



R/L 



Fig. 8. 1 Networks, transfer functions, and s-planes of passive networks 
used. The network of Fig. 8.1 (c) cannot ordinarily be used in the main 
transmission path unless there exists a parallel forward path through which 
dc components may pass. 
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TABLE 8-1 
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Definitions of the transfer functions of Fig. 8-1 for algebraic manipulation. 
The subscript d refers to poles or zeros resulting from a lead network. 

The subscript g refers to poles or zeros resulting from a lag network. 
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8.2 Investigation of Some Active Networks 



Not all active networks will yield complex poles or zeros. It is 
therefore necessary to investigate several combinations of amplifiers, 
summers, and passive networks to determine some of the configurations 
which will produce the desired complex poles and zeros. 

Two basic network configurations, one a parallel feed forward and 
the other a feedback, were selected for investigation. Solutions were 
obtained by placing a lag or a lead transfer function and amplifier in 
each path available. The overall transfer function of the network was 
then determined and investigated for the existence of complex poles or 
zeros . 

Fifteen different networks were investigated using this plan. Some 
were unsuccessful in that they did not yield the desired complex poles 
and zeros. Only those producing the desired result are treated here. 
There are, undoubtedly other simple active networks which will produce 
complex poles and zeros. 

Lead-lead feed forward difference compensator 

Fig. 8.2 is the block diagram of this compensator. The transfer 
function of Fig. 8.2 is 







(s + Pd^) + Pd2^ 



(8-4) 



For factoring by root locus, the numerator of Equation 8t4 may be 
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manipulated into the form 



( s + e^^Hs + p^,) 



(8-5) 



Note that Equation 8-5 is a zero degree root locus. The possible root 
loci arising from Equation 8-5 are Figs . 8.3, 8.4, 8.5, and 8 . 6 . It 
can be seen from Figs. 8,5 and 8.6 that complex mots of Equation 8-5 



will obtain for certain values of — — . These complex roots will become 

^2 

complex zeros in Equation 8-4 and result in Equation 8-4 having the form 



V 



V. 



(Af - A^) (s +cr+ j W ) (s +<r- j 6)) 

(^+Pdi» + 



(8-6) 



This transfer function is of further interest because p and p , the 

dl dZ 

real poles, may be placed far out in the left half plane if desired. 

There can be no complex root resulting from the root loci of Figs. 
8.5 or 8.6. Therefore they are of no further interest in this investigation. 
Lead-lead negative feedback compensator . 

Fig 8.7 IS the block diagram of this compensator. The transfer 
function of Fig. 8.7 is 




(s+ ( s + p^) 



( s + p { s + p ( s + (s + 



(8-7) 



For factoring by root locus, the denominator of Equation 8-7 may be 
manipulated into the form 
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Fig. 8.2 Block diagram of the lead-lead feed forward difference 
compensator . 



Fig. 8.3 




V 



Fig. 8.4 
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Fig. 8,5 A possible root locus of Equation 8-5. 




Fig. 8,6 A possible root locus of Equation 8-5. 
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Fig. 8.7 Block diagram of the lead-lead negative feedback compensator. 




Fig. 8.8 A possible root locus of Equation 8-8. 
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( 8 - 8 ) 



Aj ( s + (s + 



(s + p^j) (s + p^2> 



Fig. 8 8 IS the root locus of Equation 8-8. Any locations of p , p,^, 

dl 



a or will produce the locus of Fig . 8.8 unles s j [p 



dl 



d2 



d2 K rdi 



or 



> 



d2 



The complex roots which obtain for certain values of A 



will become complex poles in Equation 8-7 and result in that equation 
having the form 



V 



o 



( s + 






Pd2* 



(8-9) 



X (A^ + 1) ( s + <3" + j ( s +cr-j W) 

This transfer function is of further interest because one of the real 
zeros of Equation 8-9 may be placed close to the origin while the other 
may be placed as far out in the left half plane as desired. Note from 
Fig. 8.8 that the complex poles will always be in the left half s -plane. 
Lag-lag feed forward difference compensator . 

Fig. 8.9 is the block diagram of this compensator. The transfer 
function of Fig. 8.9 is 



V 

o 



1 =<1 ( = + (= + Pg2»- A2=<2 (s+ ZgjHs + Pg|) 



g2 



( s + Pg) ( s + 



(8 



For factoring by root locus, the numerator of Equation 8-10 may be 
manipulated into the form 
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Fig. 8. 



Fig. 8. 



Fig. 8. 1 
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9 Block diagram of the lag- lag feed forward difference 
compensator . 




10 A possible root locus of Equation 8-11. 
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Fig. 8.12 A possible root locus of Equation 8- 1 1 . 
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Fig. 8.13 A possible root locus of Equation 8- 1 1 . 
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( 8 - 11 ) 






g2' 



Note that Equation 8-11 is a zero degree root locus. The possible root 

loci arising from Equation 8-11 are Figs. 8.10, 8.11, 8.12, and 8.13. 

It can be seen from Figs. 8.10 and 8.11 that complex roots of Equation 

Aj o(i 



8-11 will obtain for certain values of 



A^C<2 



. These complex roots will 



become complex zeros in Equation 8-10 and will result in the equation 
having the form 

(Ai (s +a-+}^) (s +cr- j to) 



V 



V. 



( 8 - 12 ) 



(s+pgi)(s + pg2) 



This transfer function is of further interest because p . $ and p the 

gl g2 

real poles of Equation 8-12, may be placed fairly close to the origin in 
the left hand s-plane. 

There can be no complex roots resulting from the root loci of Figs. 
8.12 and 8.13. Therefore they are of no further interest in this investi- 
gation. 

Lag- lag negative feedback compensator. 

Fig. 8.14 is the block diagram of this compensator. The transfer 
function of Fig . 8 . 1 4 i s 



V 



Zg^)( s + 



g2' 



(8-13) 



1 (s + pgj) ( s + pg^) + AjC<jC <2 (s + Zgj)(s+Zg 2 ) 

For factoring by root locus, the denominator of Equation 8-13 may be 
manipulated into the form 
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Fig. 8.14 Block diagram of the lag-lag negative feedback compensator. 




(T 



Fig. 8.15 A possible root locus of Equation 8- 14. 
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AjC^lO^Z ( s + ( s+ 



1 



(8-14) 



(s + Pgl)( s + Pg^) 

The root locus of Equation 8-14 is Fige 8. 15. It can be seen from Fig. 
8.15 that complex roots of Equation 8-14 will obtain for certain values 
of A c( cC . These complex roots will become complex poles in Equation 

1 X tu 

8-13 and will result in that equation having the form 



V 

o 



V. 

1 



C<1 ( s + ^gj) ( s + p^^) 

1) ( s +CT+ j03) ( S +cr- joj) 



(8-15) 



This transfer function is of further interest because one of the real zeros 
of Equation 8-15 may be placed near the origin while the other may be 
as far out in the left hand s-plane as desired. The complex poles will 
always be in the left hand s-plane. 

Note that Fig. 8.15 will be the root locus of Equation 8-14 unless 
jpg^l'^ |^gl| ^ 2 I " neither of these cases could com- 

plex roots be obtained from Equation 8-14. They are therefore of no 
interest in this investigation. 

Lag-lead feed forward summing compensator . 

The block diagram of this compensator is Fig. 8. I 6 . The transfer 

function of Fig. 8. I 6 is 

A (s + 2 ) ( s + p ) + A ( s + -2 ) {s+ p ) 

^ = — i s s! ? i 8 (8-16) 

1 (s + p^) ( s + p ) 

For factoring by root locus, the numerator of Equation 8-1 6 may be 
manipulated into the form 



148 




V 

o 



Fig. 8.16 Block diagram of lag- lead feed forward summing compensator. 




Fig. 8.17 A possible root locus of Equation 8-17. 
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(8-17) 



A o(( s + Z ) (s + p ) 

I g Q 

A ( s + Z ) ( s + p ) 
d g 



Fig. 8.17 is the root locus of Equation 8 - 17 . It can be seen by inspec- 
tion of Fig. 8. 17 that complex roots will obtain for certain values of 




so long as 1 z > 
I g| 




and 



> 




These complex 



roots will become complex zeros in Equation 8-1 6 resulting in that equa- 
tion having the form 



V A o(r+ A ) { s +CT+ j M ( S + cr - j u» 

-i £ (8-18) 

i (s + p ) (s + p ) 

g d 

This transfer function is of further interest because one of the real 
poles of Equation 8-18 may be placed close to the origin in the left hand 
s-plane while the other may be placed as far out on the negative real axis 
as desired. The complex zeros will always be in the left hand s-plane. 
Lag-lead positive feedback compensator; lag in the feedback path. 

The block diagram of this compensator is Fig. 8.18. The transfer 
function of Fig. 8.18 is 



o 



( s + Z ) ( s + p ) 
a g 



(8-19) 



i (s + p ) { s + p ) - A cC( s - Z )(s + -Z ) 
g d 1 g d 

For factoring by root locus, the denominator of Equation 8-19 can be 

manipulated into the form 



A^ cC (s + Z ) ( s + Z,) 
1 g d 

( s + p ) ( s + p ) 
g d 



= 1 



( 8 - 20 ) 
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Note that Equation 8-ZO is a zero degree root locus. Possible root loci 
resulting from Equation 8-20 are Figs . 8.19, 8.20, 8.21, and 8 . 22 . It 
can be seen by inspection of Figs. 8.21 and 8.22 that complex roots of 
Equation 8-20 will obtain for certain values of cC . These complex 

roots will become complex poles in Equation 8-19 and result in that equa- 
tion having the form 



o 



( s + Z ) ( s + p ) 
d ^ g 



{ 1 - AjOC) { s +(T+ j ( s +CT - jUi) 



( 8 - 21 ) 



This transfer function is of further interest because z and p , the real 
zeros of Equation 8-21 may be placed as close to the origin of the s-plane 
as desired. There can be no complex roots resulting from the root loci 
of Figs. 8. 19 and 8,20. Therefore, they are of no further interest in 
this investigation. 

Lag- lead positive feedback compensator; lead in the feedback path. 

Fig. 8.23 IS the block diagram of this compensator. The transfer 
function of Fig. 8.23 is 



o 



( s + Z ) ( s + p ) 
g d 



( s + p ) ( s + p ) - A cC(s + Z- ) (s + Z ) 
e d 1 d g 



( 8 - 22 ) 



For factoring by root locus, the denominator of Equation 8-22 may be 
manipulated into the form 



A c(( s + Z ) { s + Z ) 

1 u g_ 

( s + p ) ( s + p ) 
d g 



(8-23) 
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Fig. 8.18 Block diagram of the lag-lead positive feedback compensator. 
Lag in the feedback path. 




Fig. 8.19 A possible root locus of Equation 8-20 . 
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Fig. 8.20 A possible root locus of Equation 8-20. 
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Fig. 8.21 A possible root locus of Equation 8-20 . 
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Fig. 8.E3 Block diagram of the lag-lead positive feedback compensator. 
Lead in the feedback path. 



154 



Equation 8-23 is identical to Equation 8-20. Therefore the root loci 
of Equation 8-23 will be those of Figs . 8.19, 8.20, 8.21, and 8.22. 
Thus, Equation 8-22 will be of the form 



V 

o 




(s + & ) ( s + p ) 

i d 

( 1 - Ap<) ( s +o~+ jW) { s +cr- jb>) 



(8-24) 



The same comments may be made for Equation 8-24 as for Equa- 
tion 8-21 except that z and p,, the real zeros of Equation 8-24 be placed 

g d 

far out in the left hand s-plane if desired. 
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8 . 3 Algebraic Equations for the Loci 



The complex zeros generated by the active compensators all 
lie on circles. For some purposes the equations of these circles 
are of interest. The derivations are simple and are not given 
here. The results fall naturally into two sets of algebraically 
similar equations, which are listed in Tables 8-2 and 8-3. 
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TABLE 8- 




TABLE 8- 



9. MANIPULATIONS AND DESIGN 



9 . I Introduction 

The equations of the root locus circles, while of definite interest, are 
not suitable to work with when designing active compensators . A convenient 
form may be obtained by manipulation of the equations of Table 8-2 and 8-3. 
For complex zeros the relationship of Table 8-2 may be manipulated into 
the form: 

(9-1) 

V^2 ‘Pl^ +2x (Zj P2-^zPl^ P 1 P 2 V ^ 1 ^ 2 ^P 2 'Pi^ " ° 



for use in design. Table 9. 1 gives the particular form of Equation 9-1 
applicable to each of the complex zero producing compensators. It will be 
shown that four of the six variables in Equation 9-1 may be chosen to meet 
the needs of a specific problem. Equation 9-1 will then reduce to an equa- 
tion in two unknowns which may be chosen so as to satisfy that equation and 
the requirement of physical realizability. 

A significant simplification of Equation 9-1 occurs if: 



P] = P2 

Zi pj i K, 



(9-2) 



Substituting Equation 9-2 into Equation 9-1 and simplifying yields 

x^ + y^ = Kj (9-3) 

Note that due to restrictions imposed by Equation 9-2, Equation 9-3 can be 
used only in those active compensators which contain two lag or two lead 
networks . 



159 






CO 

W 

O 

Oh 

< 

w 



H iz; 
o o 
o 

pc; 



H 
< 
oc; P 

o a 

H w 

CO 
Oh S 

W O 

S o 

p J 

2 



<r 

I 

<r 



(M 



0? 



(M 

TJ 

a 

4- 

co 



C«' 

4- 

03 



II 






II 


1— 1 
T> 


Ph 

4“ 


t— 1 

Ph 


?F- 


Ph 


03 


+ 


4- 


4* 

03 


(-— 1 


03 


03 


^*rs] 




(30 

CW 


"^(30' 




4- 


W ' 


03 


4“ 


4“ 


+ 




03 


03 


03 

< 


'iT 




<r 


<T 








fv] 

Pi 



a 



W) 









Pi 



Ph 

4- 

0) 

C<1 

4* 

03 



2 

o 

I— I 

H 

<! 

P 

a 

w 

o 

2 

HH 

2 

cci 

O 



Pd 

O 

< 

CO 

2 

W 

P4 

2 

o 

u 



<M 

/t3 

(X 

S 

(M 



tti 



(M 



IM 

pc, 

X 

(NJ 



(\3 

,00 



rs] 

>S 


1 

i-H 




rs] 


r— 1 

(30 


4" 


rtT 




4* 










Ph 


— ' 


(M 








o? 






1— H 




r— I 


1— 1 






Ph 

1 


P? 


TJ 

Ph 


Ph 

1 


Ph 

II 



Q? 



I 






T3 

ctJ 

0 

VM 

TJ 

(U 

(U 

T? 

nJ 

<u 

»-H 

1 

T> 

<tJ 

0) 



rs] 



II I 

fs] 

CM^ P? 



rs] 

P< .^d 



tU cU 



Ph 

4- 



W) 



tn 



rs] 

dfi 



cki 

I 



(30 



Ph 

I 



(30 



(30 

Ph 

I 

rs] 

Ph 

rs] 

(30 

D« 



U 

o 

ni 

03 

p 

(U 

Ph 

s 

o 

u 

0) 

u 

p 

<u 

u 

<u 

SM 



m;i 



TJ 

<t3 

0 
sm 

T) 

0) 

<U 

sm 

(30 

03 

I— ^ 

1 

(30 



rs] 

(30 



c« 



(30 



ti;i 



(30 



o 

4-> 

oJ 

03 

P 

<U 

Ph 

s 

o 

o 

03 

U 

P 

0) 

u 

03 

SM 



TJ 

Ph 

K 

rs] 



rs] 



rs] 



Ph 

1 

c? 



Ph 



Ph 

4* 



(30 



c« 

cm" 

II 



(30 






m 

I 

c«'^ 



(30 



TJ 


1 


cti 


(30 




c« 


Ph 




1 






Ph 


(30 


^TJ 

Ph 



T? 

OJ 

O 

Sf-I 

03 

03 

SM 

TJ 

(t3 

03 

f— I 

(30 



o 

cxJ 

03 

P 

03 

Ph 

6 

O 

u 

(30 

p 



r6TT 



TABLE 9. 



For complex poles ^ the equations of Table 8-2 may be manipulated into 
the form 

(2, + *2 - P] - Pj > ) + 2X(pj . gj Zr^) + 

*1 *2 ‘ P] Pz> • P] Pz < * *2> = “ 

for use in design. Table 9-2 gives the particular form of Equation 9-4 
applicable to each of the complex pole producing compensators. It will 
be shown that four of the six variables in Equation 9-4 may be chosen to 
meet the needs of a specific problem. Equation 9-4 will then reduce to an 
equation in two unknowns which may be chosen so as to satisfy that equa- 
tion and the requirement of physical realizability. 

A significant simplification of Equation 9-4 occurs if 

Zj z, = Pi P, ^ (9-5) 

Substituting Equation 9-5 into Equation 9-4 and simplifying yields 

- ^2 (9-6) 

Note that due to restrictions imposed by Equation 9-5, Equation 9-6 can 
be used only in those active compensators which contain a lag and a lead 
network. 
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9 . 2 Construction of the Locus of Complex Poles or Complex 
Zeros Producing a Constant Phase Angle at a Point in the 
S“plane, 

It has been shown by Carpenter (19) that if the angle at a 
point in the s-plane due to an array of poles and zeros be known, 
it is possible to construct the locus of complex conjugate poles or 
complex conjugate zeros which will cause this point to be on a 
root locus o This locus of complex poles or zeros proves to be an 
arc of a circle. Carpenter* s derivation is as follows: 

Figure 9. 1 shows the necessary construction in proving the 
locus of complex conjugate poles or zeros, R and R’, producing a 
constant phase at a point p, in the s-plane, is an arc of a circle 
passing through that point. The center of the circle lies on the 
real axis . 

From the figure: 



b = 90° +0C 

a + b = ^+ 2o<+ 180 =M=angle contributed by (9-6) 
complex pair at point p. 

But ^ + 2 cC -f ^ =180 
Subtracting 9*7 from 9.6 

360 «=■ M = ^ = a constant 

The vertex of therefore, lies on an arc pp* , of a circle 
whose center is on the real axis, by a fundamental theorem of 




plane geometry. 
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Fig. 9-1 CONSTRUCTION NECESSARY IN PROVING CARPENTER'S WORK 



Since the arc always intersects the real axis, the arc 
may be quickly constructed by first locating that intersection. 
The intersection is found by constructing the angle, 



at the point P. 

It can also be shewn that the construction angle at the point 
P which locates the center of the circular arc on the real axis is 



The necessary construction to show this is Fig. 9 -2. F is the 
intersection of the circular arc with the real axis defined by the 
angle 0 . HG is the perpendicular bisector of PF . Since the center 
of the circular arc is on the real axis, the intersection G must be 
the center of the circular arc, by plane geometry. From Fig. 9-2; 



0 = M/2 



(9-8) 



r = M 



(9-9) 



0 + g = 90° 

20 + 2g = 180° 



( 9 - 10 ) 



(9-11) 



2g = h 



( 9 - 12 ) 



Subtracting Equation 9.5 from 9.4 obtains 



20 = 180° - h 

h + r= 180° 

P = 180° - h 



(9-13) 



(9-14) 



(9-15) 



Subtracting Equation 9-15 from 9-13 
20 = r = M 



1 65 
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Fig. 9-2 CONSTRUCTION NECESSARY IN LOCATING THE CENTER OF THE CIRCULAR ARC 



f rom Equation 9 »Vj thus proving the statement. 

With the center of the circular arc and its point of inter- 
section with the real axis known from 0 and K, the arc may 
be easily constructed. 

This construction provides a sinple means of finding where 
on the s-plane complex poles or zeros may be placed in order to 
force a selected point to be on a root locus. This technique will 
not always be needed. For example when it is desired to cancel 
existing complex poles with complex zeros, the desired location 
of the complex zeros is specified by the location of the complex 
poles. However, in the more general cases, this construction 
becomes a powerful tool in problem solution. 

Note that the circular arcs resulting from this construction 
will tell the designer where complex compensating poles or zeros 
should be placed in order to cause a desired point to be on a root 
locus. The circular equations for the active compensators of 
Section 9- 1 then allows one to place complex zeros at the required 
locations . 
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9.3 A Technique for Compensation Using Active Networks 



Relocation of the root loci of the characteristic equation so that 
they pass through a specified point is not in itself any great problem. 
Simultaneously meeting the requirements of dominance, steady state 
accuracy, and transient response is, however, quite another matter. A 
considerable number of techniques exist for the simultaneous solution of 
the steady state accuracy and root location problems. See for example 
ZO, 21 . These techniques apply only to compensation with negative real 
poles and zeros and do not easily extend to the use of complex poles and 
zeros. The technique presented here requires some intuition and trial 
and error as thus far no more straight forward approach has been found 
For the trial and error phases, the Electro-measurements, Inc. , ESIAC 
algebraic Computer, Model 10 was a valuable tool. 

The techniques of solution, whether for complex poles or complex 
zeros are practically identical. Therefore they will be presented 
separately to avoid any possibility of confusion. 

(A). Solutions requiring complex zeros. 

a. Select the compensation path. Then determine the open loop 
transfer function and manipulate the characteristic equation into the 
form: 



s”,tr 



( s + Z , ) 

X 



( s + pj) 



G = - 1 
c 



(9-17) 
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b. Plot the z , p. , and s on the s-plane. 

Co Select the real poles of G so as to be most convenient. 

c 

Convenience will be determined by the particular problem. Plot these 
poles on the s-plane. This will usually determine which of the complex 
zero producing active networks is to be used. 

do Select the desired locations, P and of two complex con- 
jugate closed loop roots o Plot P and P’ on the s-plane. As mentioned 
previously, it is hoped that these roots will prove dominant, but there is 
no way of ensuring this. 

e. Calculate the phase angle at P or P' due to z^ , p^ , s^ , and 
the real poles of . Then construct the locus of complex zeros in 

the manner described in Section 9*2p which will result in P and P* being 
on a 180 degree locus. 

f. Choose two complex conjugate points, Q and Q' , on this con- 
structed locus. These will be the complex zeros of G^ . 

g. Determine the root locus gain, which will cause roots 

to be at P and P' with the z , p. , s , and the chosen real poles and 

1 J 

complex zeros of G^ . 

h. Substitute K and the real poles and complex zeros of G 

R ^ 

into the open loop transfer function to determine if steady state accuracy 

specifications are met. If not, choose two new locations of the complex 

zeros of G and repeat steps f and g. 
c 

i . . Determine x and y from the chosen locations of the complex 

zeros of G . 

c 
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j . Substitute x and y and the selected real poles of G into the 
appropriate equation of Table 9-1 to determine the algebraic relation- 
ship between the real zeros of the compensator passive networks. 

k. Construct a graph of this equation or evaluate it for several 
values of one variable and select any convenient, physically realizable 
values for the zeros of the passive networks. The z/p ratios of lead net 
works should be ^=0.1. 

l. Substitute the real poles and zeros of the passive networks into 
the numerator root locus equation of Table 9-1 corresponding to the 
active network chosen. Plot this equation and the selected complex com 
pensator zeros on another s-plane. 

m. Determine the value of A /A which will cause the roots of 

1 ^ 

the numerator of G to lie at the selected locations for the complex 

c 

zeros . 



n. Using the results of steps g and m A^ and A^ may be calcu- 
lated and the compensator design is completed. 

o. The transient response must be computed and compliance with 
specifications determined. 

(B). Solutions requiring complex poles . 

a. Select the compensation path. Then determine the open loop 
transfer function and manipulate the characteristic equation into the 
form: 

n 

K TT (s + Z. ) 

2-- G = -l (9-18) 

.ti" ( s + p.) 

1=1 J 
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b. Plot the z. , p , and on the s -plane. 

c. Select the real zeros of G so as to be most convenient. 



c 

Convenience will be determined by the particular problem. Plot these 
zeros on the s~plane. This selection will determine the particular active 
network to be used. 

do Select the desired locations, P and P’, of two complex conju- 
gate closed loop roots. Plot P and P' on the s-planeo As mentioned 
previously, it is hoped that these roots will be dominant, but there is no 
way of ensuring this. 

e. Calculate the phase angle at P or P' due to z^, p^ , s^ , and 
the real zeros of G . Then construct the locus of complex poles in the 



manner described in Section 9*2, which will result in P and P' being on 



a 180 degree locus. 

f. Choose two complex conjugate points on this constructed locus. 
These will be the complex poles of G^ . 

g. Determine the root locus gain, K , which will cause roots 

R D 

to be at P and P’ with the z„, p., s , and the chosen complex poles and 

1 J 

real zeros of G . 



c 

h. Substitute K and the complex poles and real zeros of G 

R J-/ ^ 

into the open loop transfer function of the system to determine if steady 
state accuracy specifications are met. If not, choose new locations for 
the complex poles of G and repeat steps f and g. 

i. Determine x and y from the chosen locations of the complex 



poles of G 
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j . Substitute X and y and the selected real zeros of G into the 
appropriate equation of Table 9-2 to find the algebraic relationship ber 
tween the real poles of the passive networks. 

k. Construct a graph of this equation or evaluate it for several 
values of one variable and select any convenient^ physically realizable 
values of the passive network poles. The z/p ratio for lead networks 



1. Substitute the values of the poles and zeros of the passive net- 
works into the denominator root locus equation of Table 9-2 correspond- 
ing to the active network chosen. Plot this equation and the selected 
complex compensating poles on another s-plane. 



denominator of to lie at the selected locations of the complex compen- 
sating poles. This completes the design of the compensator. 

n. The transient response must now be computed and compliance 
with specifications determined. 



should be 




m. Determine the value of which will cause the roots of the 
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9 . 4 Numerical Examples 



This section applies the techniques of solution to further clarify 
their use. The examples presented herein are straightforward so as to 
illustrate the use of active networks as compensators without becoming 
unnecessarily involved in other considerations. It is assumed that active 
compensation was necessary in these examples due to unspecified condi- 
tions . The subsections of each numerical example correspond to the steps 
in the applicable technique of solution. 

Example 1 



Given: 



The feedback control system of Fig. 9.3 =5.0. 

Requirements : 

Use cascade compensation. Desired root location to be such that 
0.5— ^ — 0.7, settling time =1.0 seconds . not to be reduced . 

Solution: 



a. The compensated open loop transfer function of Fig. 9.3 is 

K 



F = 
o 



= G 



(9-19) 



s ( s + 5s + 100) 

The characteristic equation of the compensated system is 



1 + 



K 



( s + 5s + 1 00) 



= 0 



( 9 - 20 ) 



b. The z , p , and s”^ of Equation 9.19 are plotted on Fig. 9.4. 

1 j 

c Select D = - 25, p = - 30, so as not to add any more poles 
^dl d2 

in the region near the origin, p, .. and Pj^ also plotted on Fig. 9.4. 
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Fig. 9.3 BLOCK DIAGRAM OF THE FEEDBACK CONTROL SYSTEM FOR EXAMPLE 1 . 



r 

j 




d. Select the points s= - 5 + j 7 to be P and P'^ the roots of the 
characteristic equation which will meet the requirements. 

e . F rom Fig. 9-4, the lead angle required at P is 307 . 3 degrees . 
The construction angles at P are 

0 = 153.65 degrees 
n = 307 . 3 degrees 

The arc of the circle on which complex compensating zeros can be located 
was constructed on Fig. 9-4 using these angles. 

I 

f . On Fig. 9-4, choose s = + j 8 . 8 to be and Q ^ ^ the locations 
of the complex compensator zeros. 

g. The root locus gain at P using all poles and zeros is 

h. Substituting this value of gain and the chosen compensator poles 
and zeros into the open loop transfer function obtains K =3.51. This 
does not meet requirements. Therefore, choose two new complex compen- 
sator zeros a little closer to the desired root locations. Let these new 
locations be s = -2.5 + j 8.4, and Q’ 2 ’ With these new zeros, the 
root locus gain at P is 

Kr^ =6.578 

Substituting this root locus gain and the new complex zeros into the open 
loop transfer function yields K =6.72. This meets the requirements. 

1 . From the chosen locations of the complex zeros, 

X = -2.5 
y = + 8 . 4 
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j. Svxbstituting these values of x and y, p and p into the working 

dl dZ 



equation of the lead-lead feed forward difference compensator of Table 9-1 
yields 



Z,, = 140.34 
dl 



+ 5446 

dZ 



135.34 - Z 



d2 



( 9 - 21 ) 



The algebraic relationship between z and z 

dl d2 

k. Several values of z^^ were substituted into this equation and 

z_ = -2.5, z__ = -3.0 were selected as the zeros of the lead networks . 
dl d2 

l. Substituting the chosen z , z , p , and p into the numerator 

dl dZ dl dZ 

root locus equation of the lead-lead feed forward difference compensator of 
Table 9-1 yields 



A ( s + 2.5) ( s + 30) 

{ s + 3 ) ( s +■ 25) " ^ 

Equation 9-22 and the chosen complex compensator zeros are plotted on 
Fig. 9.5. 



m. From Fig. 9.5, A /A = 0,8368 will cause the numerator zeros 

± 

to be at the chosen complex conjugate locations. 

n. The values of A and A may now be computed from the equations 

1 ^ 

A /A = 0.8368 
1 ^ 



(9-23) 



K (Aj - A^) = - 6578 

Since K is a variable gain element, any values of K, A^ , and A^, which 
satisfy Equations 9-23 will yield satisfactory results. 
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o. The characteristic equation of the compensated system may now 
be written as 



6578 (s+Z.5-fj8.4)(s4-Z.5-j8.4) 
s (s^ + 10s + 100 ) ( s + 25) { s + 30 ) 



(9-24) 



The roots of this equation were obtained by calculation on the Control Data 

Corporation l604 digital computer as 

s= - 5. 11 + j 7.26 
s= - 5. 11 - j 7.26 
s= - 4.69 + j 11.69 
s= - 4 . 69 j 11.69 
s= -40.4 

The roots as s .11 + j 7.26 and s = - 5.11 - j 7.26 correspond to those 
selected early in the solution. There is some error due to the graphical 
techniques required in the solution, however the agreement is close enough 
for engineering work. The transient response may now be calculated and 
final agreement with requirements determined. 

Epcample 2 . 



Given • 

The feedback control system of Fig. 9.6. = 100 

Requirements 

Use cascade compensation. Desired root location such that 0 . 5 
^ 0.7, settling time =1.0 seconds. not to be reduced. 

Solution: 

a. The compensated open loop transfer function of Fig. 9--6 is 



F 

o 



K 



{ s + 1) 
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(9-25) 
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The characteristic equation of the compensated system is 



1 + 



K 



s { s + 1) 



y- G = 0 

Z c 



( 9 - 26 ) 



b. The z^y p j , and of Equation 9-26 are plotted on Fig. 9.7. 

c. Select = -40, p^^ = - 50 so any root loci from these poles will 

not affect the region where the desired root will be located, p and p are 

dl dZ 

also plotted on Fig .9.7. 

d. Select the points s = -4+ jBtobeP and P' , the roots of the 
characteristic equation which will satisfy the requirements. 

e. From Fig. 9.7, the lead angle required at P is 203.8 degrees. 

The construction angles at P are 

0 =101.9 degrees 
P = 201.8 degrees 

The arc of the circle on which complex compensating zeros can be located 
was constructed on Fig. 9.7 from these angles. 

f. Choose the points s=-3.2+ j2tobeQ and Q' , the locations 
of the complex compensator zeros. 

g. The root locus gain at P and P* using all poles and zeros is 

K = 15,316 

R j-j 

h. Substituting this value of gain and the chosen poles and complex 
zeros of the compensator into the open loop transfer function obtains 



K =108.8. 

V 

This meets requirements. 



182 



i. F rom the chosen locations of the complex compensating zeros* 

X = - 3,2 
y = + 2.0 

j . Substituting these values of x and y, p and p^^ into the working 

di d2 

equation of the lead-lead feed forward difference compensator from Table 
9-1 yields 



= 169.424 
d2 



- .08405 
dl 

Z^, + 175.824 
dl 



the algebraic relationship between z and z . 

dl d2 

k. Several values of z,^ were substituted into this equation and 

dl ^ 

Z-, = -30, Z - = - 34.95 were selected as the zeros of the lead networks . 
dl d2 

l. Substituting the chosen z ,, , P P into the numerator 

dl d2 dl d2 

root locus equation of the lead-lead feed forward difference compensator 
on Table 9-1 yields 



Aj (s -h 30) (s 4 50) 

( s 4 34.95 ) ( s 4 40) 



= 1 



(9-27) 



Equation 9-27 and the chosen complex compensator zeros are plotted on 
Fig. 9.8„ 

m . From Fig. 9 . 8 A /A = 0.9314 will cause the numerator zeros 

J. 

to be at the chosen complex conjugate locations. 

n. The values of and may now be computed from the equations 



A /A = 0.9314 

J. 

K (A^ - A^) = Kri^ = 15.316 



(9-28) 
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o 




Since K is a variable gain element, any values of K, , and , which 
satisfy Equations 9-28 will yield satisfactory results. 

o. The characteristic equation of the compensated system may now 
be written as 

15.316 ( s + 3.2 + j 2) ( s + 3,2 - j 2) ^ ^ 

s ( s + 1)^ ( s + 40 ) ( s + 50) 

The roots of this equation were obtained by calculation on the Control Data 

Corporation l604 digital computer as 

s= - 3.39 + j 5.4 
s= - 3.39 - j 5.4 
s= - 4.33 
s= -19.73 
s= -61.06 

The roots as s = -3.39 + j 5.4 and s = -3.39 “j 5.4 correspond to those 
selected early in the solution. There is some error due to the graphical 
techniques required in the solution, however the agreement is close enough 
for engineering work. The transient response may now be calculated and 
final agreement with the requirements determined. 



185 



9 * 5 Summary 



For both the examples chosen, the lead-lead feed forward difference 
compensator proved to be the most practical for solution. Howeverj the 
use of the other active networks follows exactly the same procedure as 
demonstrated. 

It was found that there are some regions on the circular loci of the 
active networks at which the gain ratio A /A is extremely critical; i.e. , 

1 w 

a small variation in the gain ratio results in a large movement of the com- 
plex zeros of the active network. This sensitivity may be overcome to some 
extent by a trial and error method of varying the real poles and zeros of 
the passive networks so that the critical area of the circular locus is moved 
away from the desired location of the complex zeros. 

The active networks presented afford a simple method of generating 
complex poles and zeros. The locations of these complex poles and zeros 
may be accurately predicted and are easily varied by adjusting an ampli- 
fier gain or one or more of the passive elements of the networks. The 
varying of an amplifier gain is the more attractive of these two methods as 
it results in the complex poles or zeros moving along a circle which is pre- 
cisely defined by the passive circuit elements. 

Since the location of the complex poles and zeros is a function of 
amplifier gain a requirement for a very stable amplifier is generated. 

This requirement may be relaxed somewhat if the compensated system is 
not too sensitive to the compensator pole and zero locations. 
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The lag-lead feed forward summing compensator has characteristics 

which merit special attention. Since the K of this compensator is A + A 

c 1 Z 

the possibility presents itself of using this compensator as the main trans- 
mission path amplifier. The complex zeros resulting from the network 
could be adjusted to best suit the needs of a specific system while the over- 
all gain remained independent. 
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10. APPLICATION TO A JET AIRPLANE 



In any discussion, the argument is strengthened by an analysis 
of a practical problem, feasibility permitting. The agreement 

of calculated values and test responses of the actual vehicle 
is proof for the skeptic. This section will use the arrays presented 
in previous sections and substitute the aerodynamic stability deri- 
vatives of the North American Aviation h\c Jet Trainer (Navy 
designation T2J-1). This particular vehicle was selected as a 
model because of the wealth of flight test and aerodynamic data 

immediately available to the authors. 

The discussion will be limited to the longitudinal system 

where the aircraft, is originally in equilibrium at a level flight 
attitude. An additional comparison will show the differences 
between a two degree and three degree of freedom calculated 
estimate . 

Considering first, the following condition of flight: 



The aerodynamic stability effects were calculated and the following 
values found? 




597 

M = .6, V = ft/ sec 



o 



W = 9660 pounds, m = 300 slugs, B = 13157.7 slug-ft 



2 



Xu = - .01264 



Xw = .0384 



zu 



.0924 



Mw = - .0288 
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Mu = .0004 



Mw = - .00134 
Mq = - 1.603 



zw = - 1 . 092 
zq =- 4.00 
z = - 142.0 



M||= - 24.771 
101]^= - 1.400 
10^= - 172.00 



Xw =Xq=zw=M =10u = 10w= 10w= 10q=0 



Placing these numerical values in a determinant, the fol'owing 
form results as shown in Table 10-1. 

Evaluating the i^of this determinant, the result is ^ = ( s^+ 1 72s+ 1 . 4) 
(s^+3,5056s\l8.99122s^+.2535s+.0998) (10-1) 



This results in the following set of roots, 
s = - .01 

V elevator motion 
s= - 17K99y 

s = - .006205+ j. 07232 ( phugoid motion) 
s = - 1.7466 + j 3.9864 (short period) 



If only a short period estimate is required, then the determinant 
of Table 1C--2 is used. This assumed no speed change (u) during 
the short period and the X force equation is eliminated. 

Evaluating the A of the determinant, the results are: 
s = - .01 

(elevator response) 

s = - 171.99 

s = - 1.745 + j 3.985 (short period) 



139 



TABLE 10-1 Longitudinal System at h=30,000 ft M= 
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TABLE 10-2 Two Degree of Freedom System for T2V- 
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The short period response using the results of in Equation 10-1 



IS 



T = ■■ = 1.575 seconds per cycle 

3.9864 ^ ^ 



and the damping to 1/2 amplitude: is 



t 



• 69 



= • 395 seconds 



1/2 1.7466 



From the North American Aviation flight report NA-594-282, the 
flight test results for the T2J at conditions similar to the above 
calculated flight conditions were T= 1 . 6 seconds cycle. 



There is very good agreement between calculated and flight test 
results . 

There was insufficient information available in the North American 
report with respect to the long period longitudinal oscillation. 
Therefore, a second set of level flight condition derivativeswas 
computed using the same aircraft model. 

The conditions specified are: 



t 



1/2 



= ,4725 seconds 



h = 20 , 000 feet j, 




W =9660 pounds m=300 slugs B=13157.7 slug-ft 



2 
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For these conditions the following stability derivatives were 
found: 



XU = - .01382 
zw = - 1.425 
Mw = 8. 59 X 10 
Xw = ,0222 
zw = - 1 . 265 
Mw = - . 03365 



ZW = 0 



Mw = - .00193 



-4 



zq = - 4.60 
Mq = - 1.87 
z )| = - 130.0 
MY\ = -24.45 
10ir|= - 2.02 

10)j= - 179.0 

1 0u= 10w= 1 Owl 0q=Xw=X<^0 
M Y^= z )^= X Xy^= 0 



Thus a determinant can be formed as shown in Table 10-3. 
Evaluating the resultant A , it is found that the following roots 
exist 

s = - .0125 

V elevator motion 
s = - 178.9875 J 

s = - . 006698 , 10129 (phugoid motion) 

s= - 2.0257+’j 3.7756 ( short period motion) 

Using the results from Table 10-3, the period of phugoid 

= 62.0 seconds per cycle. 

.69 



motions is T = 2 



.10129 



The time to dampen to 1/2 amplitude is ^^12 ~ 



= 103.0 



,/2 .00670 

seconds. If a two degree of freedom estimate is desired, it can be 
determined from the above table. For the phugoid motion, the esti- 
mate can be made by considering angle of attack constant 
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TABLE 10-3 T2V-1 Longitudinal System at 20,000 feet Three Degrees of Freedom 
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TABLE 10-5 Formation of a ^ Determinant 
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TABLfE 10-4 T2V-1 Restricted to Two Degrees of Freedom at 20,000 Feet and M=.463 
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and the moment equation equal to zero in all nodes. Thus 
Table 10-4 can be made from Table 10-3o The roots from 
this determinant are: 

^=(s+178.9575) (s + .0125) (s +.0069 + j.0979) { s+ . 0069-j . 0979) 
The flight test information for the phugoid was determined 
by the authors. On 5 January 1961 ^ a T2J-1 (Navy 148204), 
piloted by McCamey and Borthwick yielded the following results 
for h=20,000 feet, M-.463; 

T -= 64 seconds per cycle 
tj ^2 ^ seconds 



These results show excellent agreement with the calculated 

estimates and prove the system to be very valid one for both long 

and short period« It reaffirms the previous knowledge that: 

lo Given reasonably accurate aerodynamic stability de.rivative, 
an accurate estimate of dynamic response can be made* 

2. For conventional subsonic aircraft, a two degree of freedom 
approximation is quite close to the three degrees of freedom 
result and should be sufficiently accurate for preliminary 
design estimates. Furthermore, factoring a second order 
polynomial is far simpler than a fourth order polynomial*, 

As mentioned in previous sections the performance function is 



stated thus 



G. K A . 

PF = J ij 

A 



(10-3) 
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It is noted that the >A. values have already been found in solving 

the natural periods of the airframe. A simple additional step 

yields the performance function desired. 

0 

For example, if the output desired the determinant 

is formed by cancelling the first row and 5th column of Table 6-1 
the results form in Table 10-5 is shown. 

The result of this determinant in algebraic form where the zero 
value terms are omitted is: 



A =s" 

15 ® 



M 






Zw -1) -Zy^ Mw 



+ s 



Xw 


Zw 


+ z 


(Mw Xu-Mw) 




- 







M y) ( Zw-Xu 






Zw-1 



(10-4) 



My^ (XwZu-XuZw) + Z (MuXw-MwXu) 

Substituting in the appropriate values for the T2J-1 at h=20,000 
feet, the following quadratic equation results 



(24.20) +s (26.900) + (.535038) 
1 5 



(10-5) 

Recalling that the direct path function for all cases in this system 



is one ( 1) , then G . = 1 . 

56 



Therefore 



© 



^56 ^5 



or 
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( 10 - 6 ) 



^ 24.20 (s+ >0203) (s+ 1.0907) 

AFe (s+.0125)(s+178.9875)(s+.006698+j 0. 1 0 1 29)( s+2 . 0257+j 3 . 7756) 

Thus the performance function for the pitch output with relation 
to elevator force is known. The block diagram of the internal 
sensing and activating equipment can be added as additional parts 
to control this multiloop system. Such an example is shown in 
Fig. 10-1 . 

The functions G1 constitutes an error sensing and signal 
amplification device such as an autopilot. G2 is the internal control 
power such as an electric or hydraulic servo mechansim G3 is the 
performance function ©/ Z\Fe . 

Assume for the moment that two input signals are desired 
in the system. For example, the aircraft may also encounter 
turbulence in the form of vertical gusts (w) as well as changes in 
elevator force ( /Vf e) . The equation for consideration of double 
inputs in general form is 

q out = . . — — G. k C ^ . q r ( 10-^) 

This last equation representing multiple inputs is shown only to 
point out the ability to use this system .for such purposes. It will 
not be used with numerical coefficients at this time, however. 



199 




FIGURE 10-1 SIMPLE BLOCK DIAGRAM OF PITCH 
REFERENCE IN AN AIRPLANE. 
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This section has demonstrated the practical use of the dis- 
cussion of previous sections. It has shown the agreement of 
calculated values and actual test. 

This particular airframe is stable, but more often supersonic 
vehicles are not. While these calculations were limited to the 
longitudinal system, the use for the lateral system involving 
the yaw and roll channel is equally valid for the appropriate 
equations of motion. The system as developed in prior sections 
shows an orderly process whereby electrical and for mechanical 
processes can be combined into one multiloop block diagram. 
From this system, with the various criteria set forth, further 
research in multiloop compensation theory can be attempted. 

The airframe problem shows one practical use for this. There 
are many others. 
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APPENDIX I 



BASIC EQUATIONS OF MOTION 



The basic concepts of mass in motion and the equations of 
motion can be derived in vector notation. This is done using mass 
particles in motion and determining their effects in two separate 
axis systems. The basic reference system shown in Fig. AI-I is 
the starred axis. This is the non- rotating reference system and 
IS also known as the Galilean Axis. The axis system is fixed in the 
rigid mass but may rotate about the fixed space axis system. It is 
the relative system. A special set of notation is used here in this 
appendix which differs from the body of this work. 

The definition of the special terms in this appendix is as follows: 
the vertorial unit along the X axis 



T - 



A 

j = the vertorial unit along the Y axis 



t ^ 



the vectorial unit along the Z axis 



i = the vectorial unit along the X axis 

A , . 

j = the vectorial unit along the Y axis 

k = the vectorial unit along the Z axis 

0 = the origin of the relative axis system 

- A 

0 = the origin of the space axis system 

the vectorial distanc 
relative axis origin 

the vectorial distanc 
mass particle in space 



A 

R = the vectorial distance from the space axis origin to the 



A 

r = the vectorial distance from the space axis origin to a 
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A ^ 

^ = the vectorial distance from the relative axis origin to 
the mass particle 

dm = the mass particle 
/O ^ 

^ = the vectorial distance from the mass center to the origi 



gin 



of the relative axis system 



A4 



rotational veloctiy vector of the relative axis system about 
the space axis; an absolute velocity vector 

^ A 

= absolute rotational velocity vector of a rotor whose axis is 
attached to XYZ axis system 

A 

L. = momentum vector for the entire mass body 

^ A 

H = moment of momentum vector with respect to the space axis 

A 

G = moment vector with respect to the relative axis 

^ A 

F = the summation of the external forces on a body 
^ = the summation of the internal forces on a body 

This appendix will derive the basic equations of motion for any 
mass moving in space. After the basic equations have been 
developed, the effects of control actions will be added The final 
equations will then be converted to standard aeronautical terminology 
which is used in the body of the work for aircraft stability considera- 
tions 

These equations have no doubt been derived several times in the 
past one hundred years. However, the author was unable to find any 
clear derivations in any of the reference material covered There- 
fore, the derivations in this appendix are submitted for the reader's 
interest and for background information 
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Consider 






ix + jy + kz 

lo = A = + ^^3 



Momentum vector for the whole body L = / dm dr 

~dt 

^ - 5 » 

= ( Vo + { OJ x^) + of rotor) 

^ A 

L = J dm dr =jdm (Vo + ( (jl>x r ) + /^ rotor) 

^ dT / 



dL 

dt 



dm d^r - + F but summation of internal forces equals zero, 
dt^ 



Hence 



dL 

dt 



dm dr = F 
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For an airframe with a large rotor (i.e. helicopter, STOL air- 
craft or GEM) consider addition of the rotor to the airframe 
using the mass center of the aircraft as the origin of the axis system. 




Vo = dR 



dt 



Consider 

d^r^ = dVo + _P— X -O-x ( x P) + P + Z 
dt^ dt 



then 
2 -?> 



d r = i Vx + j*Vy + k*Vz + i* (jfi^Vz - -^^Vy) + j* (-i^Vx- 



dt 



+ (-A- Vy - .OvVx) + 0 + 0 + 0 + 0 where /^= 0 

X ^ 



o " ^ 

Hence F = m d'^r 

i* Fx = (Vx + ^2^^ --^ 3 Vy) m 



Fy = > ( Vy + Vx - V m ? neglecting rotor effect 









k*Fz = k* (Vz Vy -jf^ Vx) m 



Therefore in scalar form 

€ 



Fx = (Vx +jf )-2 Vz Vy)m 

Fy = (Vy +,/^Vx - "^jV^) m 
Fz = ( Vz + Vy--^ Vx) m 



Concurs with Etkin 
4.4.3 

Also Duncan Page 64 
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^ ^ is fixed in the rigid body. The rotor spins at an absolute 

angular velocity of but the mass center of the rotor system is at 

point P. Therefore, and /^ = 0 . It is also assumed that 

^ c c 

_o-= constant. Therefore, the additional term-O-x ( J^x P ) 

' c 

enters the problem. 

‘ i j k ‘ 



Jl-X /^ = 



(Pi ^ 



iO-j 



Pz Pc Pc 

X ' y ’ 
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A X (B X C) - (A. C) B - (A • B) C is a basic identity, therefore 
.J^X ^ ~ ( /?) ^ 



/^C 



"r^ 

Note: coordinates i, j, k are those in the body and not the inertial 

or earth system. The values of _fl-are primed for the same 

reason. It is quite possible that will be 0 considering 

alignment of the rotor on the xz plane. 

F = Ft + F_ = m a + m_a F, = airframe F^ = rotor 
T 1 2 1 2 1 2 

In resolving this, it must be noted that the mass centers differ 
from those of the combined masses. 

It should be noted that the equations for the acceleration of a 
particle also contain rotational and relative velocity terms, however, 
f and 0 as these particles are in a rigid body. Further, the 

summation of particles about the mass center makes 

T-~^ -r> -=> ^ 

dm ( -0-x ^ ) + -N-. X { X I equal to zero . 

Moment terms may be considered using Moment of Momentum. 

From this the differential may be taken and the moment is then 

derived. Consider to be the moment of momentum of the system 

with respect to the inertial (non- rotating axis) . The derivation will 



include a rotor spinning on the rigid body. 




The absolute angular veloctiy of the body is Ul=i u| +j k LO^. 

The relative angular velocity of the body is j Ip - n( scalar) 

The angular motion of a rotor in motion, but attached to the body will be 

LAJ- -^= i ( U-'j + j (1^2 ' 2 ^ ^ ■ 

Henceforth this relative term will be noted as K (KAPPA) • 

The equation for H is: 

* - 5 > 

H = I r X dr dm where r = R +/C? 

y dT y 

Recalling that the momentum equation was L. = j dm dr 

y "dt 
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then 



= R X L + /^x dr dm 



dt 



The term under the integral will be known as 

Recalling that dr = Vo +-Q- x p + r (rotor) then 
dt 



H = 



— ^ ^ 

^ X ( Vo + (X^x P ) \ (rotor) 



dm . 



This can further be resolved into the moment of momentum about its 
own Oxyz axis;and the velocity moment with respect to the inertial 



axis . 



H = (M + m) c X Vo + P X /2-x P ) dm + ^ (K x ) dm 

-5> ^ ^ ^ 

Ho ^ X ( _fL X P ) dm 



Ho 



rel 




PP 

X ( K X ^ ) dm 






To derive moments consider recalling H = 



( r X dr ) dm 
dt 



diS' =1 

dt ) 



dr dr i dm + 
dt ^ dT“ / 



-5> 2^' 

(R + p^) X d r 

dt^ -J 



dm 



dt 



dr = -0 and recalling that F = dL> = 



dt 



dt 



2-> 

d r dm 



Then substituting these in, 
dH* =Rx dL + G where G = 



dt 



dt 






'x d^r^ 



G dt J 



dm 



However, we also noted that 



H* = R X L + 



7^ \ 

p X dr Idm 
dt 



211 



Therefore* 



^ X L + R X + dH 

dt dt dt dt 

Equating the two results of dH^' it is seen that 
. ^ dt / 

R X y^dL + G = dR x L + R\jc ^dL, + dH 
^ dt dt dt dt 

Recall that dR = Vo 
"dt 



and that 

-z5> ^ ^ 

L = (M + m) (Vo +Xlx rotor) 

/ c ■' r 



Now_^= £_ 
dt dt 



dt 



combining 



't ^ ^ 

(M + m) X Vo t Ho t Ho , 

/ c rel 



(m + M) 




X Vo + 


1? 

P X £Vo 


+_^ 


r ^ " 

Ho + Ho 

rel 




dt 




^ dt 


dt 


— 



G = Vo X (Vo +/\x ) M + Vo X (Vo +-0-x rotor) m 

^ c ^ ' c ^ c 

•’ ' — 5> 

+ ( P ^ X Vo + ( ,XLx ^ xVo) vc\ \ { p xVo+( jX-x p) X Vo) M 



+ (M + m) / X dVo + d 



Ho + Ho 



rel 



A further assumption is stated at this point. This is that the point 
of rotation of the rotor is fixed on the body and remains so. The shaft 
does not have nutation about some axis but remains fixed on a mount 
o£ the rigid body. Therefore, the ^ of the rotating body is zero. 
Also, logically, in a rigid body the mass does not shift hence 



is also zero . 
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In crossing a vector vs/^ith itself Vo x Vo the result is zero. 



Vo X Vo = 



1 J 

V V 

1 2 



k 

V, 



V V V 
1 2 3 



=t(V^V3- 

+ k (V^ V^ - V^v^ ) = 0 



Vo X Vo = 0 

Furthermore, crossing two dissimilar vectors in opposite order yields 
opposite signs and equal magnitude 

j k 

--‘<''2 /I- /I V3)^(/fv3-/5v,) 

(V, ^ V 3 ) 



Vo X /^ = 

• C 



V, 



2 3 



P\ Pi 



-5> 



/fxV?=f" ( ^3 V3 - ^3V3 ) ,f( -^,V3) +M 



Thus Vo X ^ ^ X Vo = 0 






Incorporating these two facts, the equation for 6 can be simplified to 
the following: 









+ (M + m) ( /P X dVo ) 



G = d Ho + Ho 

dT L 

-5> 

dHo = Ho + J L X. Ho in a general equation, 
dt 



21 s 






Thus it is necessary to solve dHo to find G 

^ ^ -s> 



X ( to X 1^5 ) 



Considering Ho = 
absolute velocity. 

Recalling an identity of vector algebra: 
A X (B X C) = (A* C) B - (A . B) C then 
/^x ( X Z’) = (/^./>) ( /^ uJ ) /^ 



which accounts for the 



Ho = 



2 2 2 

(x + y + z ) (i tOj + j 0^2 + tAj) - (x UJj + Y (/)2 + ^ 



^ -p 

(ix + jy + kz) 



dm 



i 



/ 2 2 , , 

(y + z ) xy - xz (Jd^ 



+ jjf(x + z ) u)2 - yzU^ - xy 

■^r 2 2 n 

+ k (x + y ) {Js^ - xz uJj yzt<^ 



dm 



The above terms can be recognized in part as moment and product 
of inertia terms. 

A / , 2 2 . 



Ix — f 

T ^ /^/ 2 2, ^ 

ly — (x + z ) dm 



Iz 



^ T/ 2 2 , ^ 

= I (x + y ) dm 






Ixy = J (xy) dm 
lyz ^ J (yz) dm 
Ixz = / (xz) dm 



2 K: 



Therefore: 



Ho 



•> — 

= i ( Ix tOj - IxyCO^ - Ixz UJ^) + j (-IxyU| + ly - lyz UJ^) 

+ ^ (Iz (jOj - Ixz LC*^ - lyz 



Considering the effects of the aircraft as a rigid body first 



— ?> -r-^> -3> ^ ^ X 

dHo = Ho +.Xi-x Ho and for a rigid body tO = -W- 
dt 

dHo = i (Ix - Ixy - Ixz ) + j (- Ixy-^ + ly-^^ -lyz 
dt 12 3 12 3 



— T. « « 

+ k (-Ixy/1- - lyz + Iz + 

X ^ J 



. 0-1 



J 

n_2 



k 

-0-3 



(Ix-Ixy-Ixz) (ly-Ixy-Iyz) ( -Ixz-Iyz+Iz) 



For any ordinary aircraft or missile, the produce of inertia terms 
Ixy and lyz can be considered zero because of the mirror symmetry 
aspect. Ixz may be zero if the principal axis is selected. Thus with 
no loss in generality, consider Ixy = Ixz = lyz = 0 

d^ [ix + Iz -O-^ -^3 - ly -^2 3 1 

dt r- » 

+ j + IxJXj^ - IzD^ J 

+ k j^iz/L^ + ly-H-j -Ol^ - Ixn^ Ci^ 
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However, including the effects of Ixz, the product of inertia term, 
which is often not zero: 

dHo = i r Ix SL - Ixz -H. Ixz /I + lzJ\^ - ly ^ 

— L 1 ^ 1 ^ ^ Z i 



4*^ Fiz -^2 “ Ixz /i-j + Iy^2*^2 " ^y^ 2 Z ^ ^2 3 



Using standard aerodynamic terminology, the following items are 
defined: 






Iz 



p=n^ 

Q -^2 

R --0-3 



Thus 




where system is chosen such that 




= 0 



-> 1 


r • * 1 


G=i 


AP - ER - EPQ + QR (C-B) I 


1 

— ^ I 


r * V 2 2 1 




BQ + PR (A-C) + E (P - R ) 




L -* 


1 

1 


|CR - EP + PQ (B-C) + EQR 



These equations concur with Duncan, Page 66, and Etkin, Chapter 4. 
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To consider the rotor effects it should be remembered that 






> ^ 

= — 0-.+ K. 'where jC is the relative motion of a rotor on the body. 



Thus for a rotor there will be a combination of terms considering 
the absolute angular velocity of the OXYZ system plus the angular 
velocity of the rotor with respect to the body. 

Thus 

dHo = Ho + t/0 X Ho whe r e{j) ^ J1-+ rC. 
dt 



It should be noted, that for a symmetrical rotor or disk the 

product of inertia terms are zero where the disk is at least a 

three bladed propeller. The two bladed propeller has a product 

of inertia dependent on its angular position. This case will not be 

considered, but the three or more bladed propeller situation will 

be used in all considerations. Hence with a dynamically and 

statically balanced rotor ^ = 0 and the product of inertia is zero. 

The moment of inertia of the disk must be considered from 

the axis of the rigid body. Therefore, the moment of inertia 

2 

must include the md consideration. 
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The moment of inertia terms for the rotors will be thus defined 

in the most general sense 
^ ^ 2 

a = Ix rotor = Ix + md^ and aligned with the reference axis 

b = ly rotor = ly + md^^ and aligned with the reference axis 

c ^ Iz rotor = Iz + md^^ and aligned with the reference axis 



G = dH 
dt 



- dHo + dH 



dt 



rel 



dt 



Recalling the derivation of moments for the rigid body without the 
rotor, the same method applies. 



218 



for a rotor only 



dHo = i a + jb ~CL^ + kc XL^ + Ho 

- 1 o j 

dt 



V ~5> -?> 

+ i a V( j + jb + kc X Ho 



Combining this to a system consisting of a rigid body plus a rotor, 
the following results occur 



>. -rs- -> 
G = dH = i 



dt 



Ix - Ixz + Ixz _T^ +1^ -O-^ ■ ly "^2 "^3 



ly 0-^^ + IxzOj^ -Iz _r^ 

— M> • 

Iz/1 - Ixz O- + ly -^1 Ix^ I.XZ /I -1^. 

3 l lz IZ Z, 



+ j 

+ k 

,y ,'y — > / 

+ iaj(^ + jbl <2 + kc l<^^ + 



i 



-5> 

j 



ill ilZ il3 



k. »K, 4 



Combining and using aerodynamic notation. 



G = r I^AP - ER + a 1^ - EPQ + (C-B) QR + ck^Q - bX^R ^ 

+ ^ |bQ + b PR + E(P^ - R^) + a ICjR - c K^p] 

+~k I^CR - EP + c + (B-A) PQ + EQR + bl^ P - aKjQ 

This result concurs with Etkin, page Il6. 

Therefore, using a standard aerodynamic notation the moment: 
G = i L + j M + kN 



and the scalar quantities are as follows 
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L = AP + a K - E (PQ +R) + QR (C-B) + c K,Q - b K.R 
1 ,32 

M + BQ + + PR ( A-C) + E (P^ - R^) + R - c\(^R 

N = CR + c + E (QR - P) + PQ (B-A) + h Q 

Thus, the force and moment equations derived in this appendix 
are applicable to a rigid body of constant mass where the body is 
in motion. It also considers the effects of 21 three bladed or more 
rotor system which is attached on some fixed point of the rigid body 
The rotor effects may be included or discarded depending on the 
magnitude of rotor inertia and the relative velocity of the rotor. 

It should be noted that the derivations are done in a general 
form, but in the concluding equations, the usual aerodynamic nota- 
tion is substituted to bring the equations to this specific problem. 
Furthermore, missile or aircraft shape and density has been 
considered symmetrical about the xz plane. Thus, the produce of 
inertia terms Ixy and lyz have been considered zero. No generality 
is lost by this and if, for some reason, the body is sufficiently 
unsymmetrical , the terms can be quickly added in. 



The following is taken from Chapter 4 of Etkin (4) . The concept 
of motion on the control surfaces is shown in the following 
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The LaGrange equation of motion in a moving reference frame is 



dt 






^ T _ V 



Zl,. 



T = kinetic energy of the system relative to the chosen frame of 
reference 



= generalized force 



W = work done on the system by the external forces which act upon 
it 

qj^= generalized coordinate 

H = aerodynamic moments 

(^ = rotation of control surface in radians 

• 2 

Kinetic energy T = 1/2 and ^ T = 

^ '■K, 



/<!K. 



For elevator motion 



= He/ (/e) + Fe/(/e) + /wi 



Fe ^ generalized elevator control force 

( ^ 

^ Wi = work done by inertia forces 

Work done by the control forces is 

Fecf(cTe) = PcT sp + J(f sj 

Fe = P dsp + J dsJ 

dc/ e d/e 

dsp 

y = gear rate of pilot’s stick to control surface 

d^ 

— = ’’gear ratio” or advantage of hydraulic cylinder to contrd surface 



Evaluating Wi 

2 

consider dFi = dr dm 
dt 

recalling from previous derivation that 

d r = dVo + ,n.x ^ ii x/^+ / + 2d%c/^ 

dt dt 



Then 



dFxi = - 



dVo + 2 Qz - 2 Ry - X (Q^ + R^) + y (PQ - R)+z(PR-Q) 



dt 

dFyi and dF z are similar in form. 

Looking at an elevator assembly, it can be seen that assuming a 
lamina in the xy-plane, the displacement in the C direction is 
amount ^ e) . 

The work done by the inertia forces is ^ ( cT e) dFzi 



dm 



cf Wi 



= I dVoz y^e dm - (PR-Q) ( ^ - RQ + P ^ y/%m 

y dt ^ 



The coriolis terms 2 Py and -2 Qx are zero as there is no work done 
in those directions. Furthermore z = 0 because of the assumed lamina. 



acz = dVoz and ac: 
“dt 




e dm = m 1 a 
e e c 



(PR - Q) j X /^e dm = P^ (PR -Q) 



A 



P product of inertia of the elevator 

X 
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cT 






By symmetry / yT y dm = 0 



J'CT e) 



Wi = - m 1 a - P (PR-Q) 
— e e c ex 

z 



I oT = m 1 a - P (PR-Q) +H + F 

e e eec e e e 

z X 

7 

I (V + m 1 a + P (PR -Q) =H + F 



eec ex 

z 



e e 



(Etkin 4 



Similarly for the rudder 



I 6 -ami - (R + PQ) P - (RQ - P) P = H + F 
r ^ r cy r r rx rz r c 



and for the aileron 



I ^ + 2 P (RQ + P) = 2 H + F 
a ^ ay a a 



8, 13) 
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APPENDIX II 



ROLL, YAW AND PITCH 

The derivation of the yaw, roll and pitch rates is given here to 
show the actual motion taking place Frequently the motion is 
considered small and small angle approximations are used. 

These will be noted at the conclusion. 

First, consider the Eulerian Axis used in the aeronautical 
field shown below. The rotation takes place in the order of yaw 
( (y ) , pitch (© ) , and roll (0 ) . 




Figure AH- 1 EULERIAN AXIS SYSTEM SHOWING ROTATION 
ABOUT THE Z, Y, AND X AXIS 
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The initial conditions at X , Y , and Z and the relations of 

o o o 

^ • • 

and 0 can be related to P, Q and R of the body axis at 
X, y, z. In the above diagram consider the angular velocity 
P which is about the OX axis. It must be the sum of the velocities 
about the and axis respectively. 



TT 



OZ^ makes an angle ^ 

TT 



+ ■© with OX. The rate thus 



cos ( — + 0) = - tpsin ©. OY is perpendicular with OX 

iLf J. 

• TT 

so © cos — - — = 0 and OX is parallel with OX or the motion 

« 0 

0 cos © = 0 . 

Thus the summation P = 0 - 4^ sin 0. 

To evaluate Q the angle YOZ must be determined. Specifically 

o 

• 

the value cos YOZ is desired. 

' o 

/\ /A 

cos YOZ = cos YOZ cos Z OZ 
o 2 2 0 

TT 

= cos ( r — - 0) COS © 



= sin 0 cos © 



Thus the contribution of OZ axis motion is ^ sin 0 cos ©. 

o 

The angle between OY and OY^^ is 0 so the motion here results as 
«> 

© cos 0. OY is perpendicular to angular motion 

0 cos ^ = 0. 

4 

Thus Q = 0 cos 0 + sin 0 cos © 

In the third angular velocity motion R the yaw motion appears as 
cos 0 cos © while the motion of pitch is © cos ( — — + 0) = -© sin 
0. The rolling term vanishes 0 cos = 0. Thus R = - © sin 0 



+ LJ^cos 0 cos © . 
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Hence the relationship of the body yaw, pitch and roll rates 
(P, Q, and R) to the designated reference system angular velocities 
Cp , © and 0 is: 

* » 

P = 0 - ip sine 
* 

Q = © cos 0 + (psin 0 cos © 

R = cos 0 cos © - © sin 0 

These results concur with Duncan page 81 and Etkin page ll6. 

It should be noted that for small angles the approximation below 
IS often used: 

0 

Q = e 



A second consideration which must be made regarding movement 
about the reference axis is the velocity or displacement with respect 
to time. It is desired to obtain the velocity with respect to the 
fixed frame or inertial axis x, j*y, and 
Consider first that in scalar values 



dx» =U 
dt 



dy^*< =V 
dt 



1 



1 



dz^ =W 
dt 



1 



It should be noted, however, that in the three position system shown 
on the previous page that 
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= V cos 0 - W sin 0 
= V sin 0 + W cos 0 
and W = - U sin *© + W cos © 





V 

3 










«2 = 




cos 


e + 


^3 


sin © 


and U^ = 




cos 


1 ^- 




sin Ip 


II 

r-H 

> 


«2 


sin 




''z 


cos jt 


Wj- 


W2 











Combining the 1 , 2 and 3 position terms 

dx'*" = U cos © cos V (sin 0 sin © cos cos 0 sin Cp) 

d^ 

+ W (cos 0 sin © cos (p-\- sin 0 sin Ip ) 
dy'l' = U cos© sin ^-1- V (sin0 sin© sin cos 0 cos ^ 

dT“ 

+ W (cos 0 sin © sin (^- sin 0 cos (^ ) 

dz-»- = -U sin © -h V sin 0 cos © -I- W cos 0 cos © 
dt 

Thus, the velocity or rate of change with respect to the inertial 
or earth’s reference axis in this case is expressed in these terms 
as derived. These velocity terms are derived in Etkin Page 102 
and conincide with the derivation shown. 

Therefore, this appendix accounts for the changes in angular and 
directional displacement of a rigid body with respect to a fixed 
reference axis system. 

22 ? 



APPENDIX III 



NOTES ON DETERMINANT 
MANIPULATION 



The body of this project contains several determinant expan- 
sions where the order is doubled. These manipulations do not 
alter the values of the determinant and the resulting values of the 
performance functions are the same. 

The simple determinant 



A1 


B1 


Cl 


A2 


B2 


C2 


A3 


B3 


C3 



evaluated will yield the following result 




A1B2C3 + B1C2A3+ C1A2B3 




-C1B2A3 - A1C2B3 - B1A2C3 

This expansion by the equations -1 + 1 = 
determinant. 



A3-1 

0 will result in the following 



A1 


0 


B1 


0 


Cl 


0 


-1 


1 


0 


0 


0 


0 


0 


A2 


B2 


0 


C2 


0 


0 


0 


-1 


1 


0 


0 


0 


A3 


0 


B3 


C3 


0 


0 


0 


0 


0 


-1 


1 



The resulting evaluation will yield the same value as result (A3-1). 
Furthermore, if the terms below the main diagonal are allowed to 
remain in their original columns in lieu of shifting one place to the 
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right as is done in the main body, the resulting ^.is still the same. 



The results of evaluating cofactors for particular performance 

functions remain unchanged. For example it is desired to evaluate 

the cofactor of the original determinant. Thus: 

1 



A = 

12 



-A1 B|i — -ei- 



A2 


B|2 


C2 


= (-1) ^2C3 - A3C^ 


A3-2 


A3 


b'3 
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In the expanded determinant the equivalent cofactor 
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The result of Equation A3- 3 equals Equation A3=2. With the 
values below the main diagonal left in their original columns, 
the following occurs: 
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A1 0 

-1 1 

A2 0 

0 0 

A3 0 

0 0 



B1 0 Cl 0 

0 0 0 0 

B2 0 C2 0 

-110 0 

B3 0 C3 0 

0 0-11 



By inspection of columns two and six, two reductions are made 
to the following; 



A 

13 



A2 0 C2 

0 1 0 

A3 0 C3 




(- 1 ) 




A3-4 



The resulting equation (A3-4) is the same as equation (A3-2). 

A similar examination using a four by four determinant 
system was made. This was done to insure that the simplified 
effects of a three by three or lower order system did not cover 
up an error in a four by four or higher order determinant. The 
same manipulations can be performed and the resultant value of 
the determinant is not altered. 
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Thus, the individual has the option of placing his values 
in either column« The location of the values is dependent on 
desired location of the pick-off point on the block diagram. 
The final result is not altered 
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